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Introduction 


Near-infrared  light  offers  unique  possibihties  for  revealing  physiological  information  from 
within  biological  tissues  without  injections,  surgery,  dangerous  radiation,  or  great  expense. 
Although  it  has  been  known  for  some  time  that  non-ionizing  near-infrared  hght  penetrates  a  few 
centimeters  into  biological  tissues,  only  recently  have  scientists  been  able  to  take  advantage  of 
this  penetration  by  taking  into  account  tissue  scattering.  Near-infrared  hght  provides  a  way  to 
quantify  the  concentration  and  the  oxygenation  state  of  hemoglobin  and  also  the  scattering 
properties  of  the  tissue  and  thus  provide  criteria  to  characterize  breast  lesions  in  vivo  without 
surgery.  This  optical  information  wiU  provide  clinicians  with  entirely  new  criteria  for  deciding  if 
lesions  are  benign  or  mahgnant,  which  may  help  reduce  some  of  the  ambiguities  and  false- 
negatives  currently  encountered  in  conventional  breast  examinations. 

OpticaUy  monitoring  physiological  changes  in  breast  tissue  also  provides  a  new  way  to  look  at 
both  disease  progression  and  prevention.  Using  near-infrared  hght  from  an  instrument  we  have 
developed,  we  have  observed  changes  in  hemodynamics  and  structure  that  are  characteristic  of 
tumors.  In  addition,  we  have  also  observed  changes  in  breast  tissue  physiology  resulting  from 
hormone  replacement  therapy,  menopausal  status,  and  estrous  cycle  stage.  Because  our 
technique  is  safe,  non-invasive  and  comparable  in  cost  to  ultrasound,  patients  can  be  studied  as 
often  as  necessary  without  risk.  Currently,  we  are  in  the  process  of  performing  clinical 
measurements  on  women  with  breast  tumors  and  on  normal  volunteers  receiving  hormone 
replacement  therapy  with  selective  estrogen  receptor  modulator  drugs  (such  as  tamoxifen). 
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Items  in  the  Original  Specific  Aims 


For  convenience,  I  will  list  below  summaries  of  the  specific  aims  of  the  original  proposal, 

“Measurement  of  Breast  Optical  Properties”: 

(1)  Optimize  the  instrumentation  and  theory  currently  used  to  extract  the  absolute  optical 
absorption  and  reduced  scattering  parameters. 

(2)  Perform  clinical  measurements  of  breast  tissue  optical  properties  in  order  to  determine 
exact  absorption  and  scattering  factors  in  different  tissue  regions. 

(3)  Compare  optical  property  measurements  to  gross  pathology,  histopathology,  ultrasound, 
and  mammography  results  in  order  to  determine  absorption  and  scattering  factors 
characteristic  of  tumor  and  normal  tissue. 


Summary  of  Progress 

Important  leaps  in  my  training  have  led  me  to  make  significant  progress  on  two  of  the  three 
specific  aims.  Enough  testing  and  optimization  have  gone  into  place  that  has  enabled  me  to 
make  significant  improvements  to  the  instrumentation  such  that  a  second  instrument  will  be 
ready  this  year  for  permanent  clinical  use.  I  have  also  performed  over  40  clinical  measurements 
that  have  led  to  two  publications  and  two  more  submitted  publications.  These  measurements  are 
of  vital  importance  since  the  provide  information  as  to  the  scattering  and  absorption  contrast  of 
near-infrared  optical  methods  in  the  breast. 


Results  of  Clinical  Measurements 
Clinical  Summary:  Normals 

I  have  been  able  to  quantify  via  absorption  the  concentrations  of  blood,  water,  and  fat  in  the 
breasts  of  over  35  healthy  volunteers.  The  addition  of  fat  has  provided  some  unique  insights  into 
breast  physiology  that  were  not  available  to  us  only  a  year  ago.  In  addition,  I  now  have  evidenee 
that  suggests  that  the  scattering  properties  of  the  breast  are  highly  influenced  by  the  amounts  of 
fat  and  connective  tissue,  and  thus  provide  complementary  information  about  breast 
composition.  This  newly  introduced  spectroscopic  information  has  provided  a  more  complete 
picture  of  the  breast,  the  first  of  its  kind. 

Current  data  have  demonstrated  several  things: 

(1)  Increased  blood  and  water  in  pre-menopausal  versus  post-menopausal  breasts. 

(2)  General  increase  in  hpids  in  the  post-menopausal  breast. 

(3)  Increased  scattering  in  the  pre-menopausal  breast,  which  is  likely  due  to  increased 
connective  tissue  density  and  decreased  amounts  of  fat. 
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(4)  Tissue  saturation  measurements  of  pre-menopausal  women  are  statistically  lower 
than  those  of  post-menopausal  women.  In  addition,  post-menopausal  women  on 
hormone  replacement  therapy  (HRT)  show  tissue  saturation  values  closer  to  the 
pre-menopausal  breast  than  the  post-menopausal  breast.  (This  saturation  is  a 
measure  of  the  hemoglobin  saturation  in  the  tissue  with  both  arterial  and  venous 
contributions). 


AH  of  the  findings  are  not  surprising,  but  until  now,  none  of  these  compositional  properties  of  the 

tissue  could  be  measured  non-invasively  in  a  safe  manner. 

This  information  is  important  for  several  reasons: 

(1)  These  are  the  first  data  anywhere  that  help  establish  a  baseline  of  “normal”  optical 
properties  for  use  in  cancer  detection  and  diagnosis.  Such  data  will  help  by  providing 
background  optical  properties  to  investigators  perusing  optical  mammography. 

(2)  Knowledge  of  “normal”  optical  properties  also  provides  potential  for  monitoring  any 
treatments  (i.e.,  HRT,  neo-adjuvant)  that  will  affect  the  composition  of  breast  tissue.  In 
addition,  we  have  demonstrated  metabolic  changes  in  the  breast  that  are  due  to  both 
endogenous  and  exogenous  hormonal  influence  (i.e.  estrogen  and  progesterone). 

(3)  For  the  first  time  we  might  possibly  quantify  the  origins  of  radiological  mammographic 
density,  which  is  thought  to  correlate  with  breast  cancer  risk.  To  date,  the  physiology 
behind  mammographic  density  is  speculated.  This  technique  could  provide  the  answers. 

These  findings  have  become  the  subject  of  two  research  papers  which  have  been  recently  been 

submitted  for  publication. 


Clinical  Summary:  Cancer 

I  also  managed  to  perform  an  interesting  measurement  that  showed  significant  changes  in  the 
breasts  of  a  patient  receiving  chemotherapy:  cytoxan  (20  mg/mL)  and  doxorubiein  (2  mg/mL). 
Prior  to  therapy,  the  tissue  saturation  of  the  tumor  was  very  low,  about  44%.  (Errors  are  difficult 
to  estimate  for  this  subject,  but  certainly  an  absolute  error  of  8%  is  a  conservative  estimate). 

Although  44%  is  somewhat  low,  low  tissue  saturation  is  representative  of  high  metabolic 
demand,  particularly  that  of  tumors.  After  the  first  bolus  of  chemotherapy  agents,  the  saturation 
ballooned  to  73%,  a  highly  significant  change.  This  large  change  in  tissue  saturation  is  evidence 
that  the  cells  are  relaxing  their  demand  for  oxygen,  implying  that  these  cells  are  dying.  Normal 
saturation  values  for  the  breast  fall  into  the  70-80%  range.  A  measurement  3  weeks  later 
revealed  that  the  tissue  saturation  remained  at  72%.  Interestingly,  the  normal  breast  showed  a 
61%  saturation  before  treatment  and  a  75%  after.  Another  measurement  3  weeks  later  showed 
that  the  tissue  saturations  of  both  breasts  were  equal  at  75%,  which  is  a  perfectly  normal  value. 
It  is  important  to  state  that  the  term  “perfectly  normal  value”  would  not  be  possible  without  the 
data  I  have  obtained  in  our  study  of  normal  subjects. 
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These  measurements  are  very  preliminary  and  of  great  value  to  the  community  of  therapeutic 
intervention.  Our  technique  monitored  the  changes  in  tumor  oxygen  metabolism.  Given  that  the 
effect  of  metabohc  decrease  was  very  high  after  the  initial  treatment  and  then  leveled  off 
indicates  that  subsequent  treatments  may  have  been  ineffective.  Perhaps  we  have  witnessed 
resistance  of  the  tumor  to  the  chemotherapy.  These  questions  are  of  critical  importance  in 
battling  cancer  and  modulating  treatment  that  is  specific  for  individual  tumor  response. 


Instrument  Improvements 

Two  technical  improvements  promise  to  forward  the  achievement  of  the  initial  specific  aims. 
We  have  incorporated  a  steady-state  spectrometer  into  the  FDPM  instrument.  The  apphcation  of 
a  simple  model  then  transforms  a  seven-wavelength  instrument  to  a  hundreds  of  wavelengths 
instrument.  The  topic  of  introducing  the  spectrometer  into  the  system  was  the  subject  of  a 
recently  accepted  paper.  Adding  the  spectrometer  should  provide  previously  unimaginable 
precision  in  determining  the  physiological  makeup  of  tissues. 

I  am  also  on  the  verge  of  completing  a  second  FDPM  instrument  that  wiU  reside  permanently  in 
the  Breast  Center  at  the  Chao  Comprehensive  Cancer  Center.  This  instrument  represents  some 
significant  technical  advances,  including  a  completely  fiber-based  probe  and  faster  measurement 
times.  Completion  of  the  instrument  is  expected  in  the  fourth  quarter  of  2000.  Personally,  I  am 
very  excited  by  this  instrument  since  it  should  become  a  means  fi'om  which  to  collect  vast 
amounts  of  data,  most  particularly  on  cancer  patients. 


Protocol  Design,  Construction,  and  Implementation 

I  have  continued  to  manage  two  research  protocols  here  at  UCI  (#95-563  &  #99-2183).  I  have 
added  some  modifications  to  these  protocols  that  should  help  us  gain  further  insight  into  the 
optics  of  breast  tissue.  One  example  is  a  measurement  we  have  started  that  will  take  place  over 
the  entire  menstrual  cycle  of  a  27  year  old  volunteer.  No  difficulties  or  comphcations  have 
occurred  during  the  courses  of  these  trials. 


Clinician  Interaction 

I  have  continued  to  work  with  clinicians  here  at  UCI.  A  recent  addition  to  the  research  team  is 
Dr.  David  Hsiang,  a  clinical  breast  cancer  surgeon.  He  and  Dr.  Butler  will  be  instrumental  in 
orchestrating  breast  cancer  measurements  at  the  Chao  Comprehensive  Cancer  Center  in  the 
upcoming  year.  I  have  had  the  opportunity  to  see  breast  cancer  in  a  more  clinical  hght  after 
consultations  with  these  physicians.  Their  interest  in  research  has  proven  invaluable  to  my 
development  as  a  breast  cancer  researcher. 


Working  with  Students 

I  am  also  continuing  to  work  with  Dr.  Tromberg’s  students.  Dr.  Tromberg  now  has  four  students 
working  in  FDPM,  and  I  am  responsible  for  training  and  guiding  them  in  their  collective  study  of 
biomedical  optics.  Ryan  Banning  and  I  have  begun  looking  at  optical  measurements  of  human 


7 


tumors  in  mice  in  order  to  better  understand  this  key  link  between  the  microscopic  tumor 
characteristics  (i.e.,  vessel  counts,  collagen  disruption)  and  the  macroscopic  optical  properties 
that  we  measure  in  vivo.  Natasha  Shah  and  Dorota  Jakubowski  and  I  have  continued  to  look  at 
optical  measurements  of  both  healthy  and  diseased  breast.  Finally,  David  Cuccia  and  I  have 
begun  some  work  that  will  correlate  optical  signals  in  tumors  with  high-resolution  MRl  images. 


Proliferation  of  Results 

Over  the  last  year  1  have  presented  my  results  at  several  meetings  and  have  published/submitted 
four  manuscripts.  These  are  listed  in  the  reportable  outcomes  section.  In  addition  to  these,  1  will 
be  attending  the  SPIE  Optics  meeting  In  January  (San  Jose,  CA). 


Attendance  at  Classes 

1  have  also  begun  to  observe  a  class  at  UCl  entitled  “Clinical  Cancer  for  Basic  Scientists.”  It  is  a 
graduate-level  course  that  introduces  cancer  researchers  to  how  clinical  practice  is  performed. 


Targets  for  Next  Year 

There  are  stiU  many  more  normals  to  observe.  Just  now,  we  are  entering  into  a  position  where 
we  can  seriously  address  issues  involved  with  cancer  diagnosis  and  treatment.  For  example,  we 
are  beginning  to  look  into  correlations  between  our  compositional  findings  and  breast  cancer 
risk.  For  this  part  of  the  study,  we  will  require  many  normal  as  well  as  cancer  patients.  .  The 
results  of  our  neo-adjuvant  study  are  highly  preliminary,  but  warrant  more  careful  study. 
Continuation  of  this  protocol  allows  us  to  verify  our  findings,  especially  with  a  second 
instrument  at  UCIMC.  We  therefore  request  another  year  of  study  to  help  mature  the  FDPM 
technique  into  a  useful  clinical  tool. 
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Appendix 


Research-Related  Accomplishments 

•  successful  measurements  on  40  subjects  (mostly  normals) 

•  shown  that  compositional  and  functional  differences  can  be  observed  in  the  breast  using 
non-invasive  optical  techniques 

•  demonstrated  a  preliminary  link  between  a  non-invasive  measurement  of  metabolism 
(oxygen  saturation)  and  hormonal  influence 

•  performed  one  of  the  first  measurements  on  the  monitoring  of  the  efficacy  of 
chemotherapy  agents  in  vivo 

•  second  generation  instrument  completion  expected  by  end  of  calendar  year 


Reportable  Outcomes:  Conferences 

[1]  A.  E.  Cerussi.  A.  J.  Berger,  F.  Bevilacqua,  D.  Jakubowski,  N.  Shah,  and  B.  J.  Tromberg,  “Non-invasive 
quantification  of  breast  tissue  physiology;  Optical  measurements  of  blood,  water,  and  lipids,”  Biomedical 
Engineering  Society  2000  Meeting,  Seattle,  WA,  Presentation  14794,  (October  2000). 

[2]  A.  E.  Cerussi.  A.  J.  Berger,  F.  Bevilacqua,  N.  Shah,  D.  Jakubowski,  J.  Butler,  R.  F.  Holcombe,  and  B.  J. 
Tromberg,  “Quantitative  non-invasive  optical  spectroscopy  of  breast  tissue  physiology,”  Optical  Imaging 
Workshop:  From  Bench  to  Bedside  at  the  NIH,  Bethesda,  MD,  P-23,  (September  2000). 

[3]  A.  E.  Cerussi.  N.  Shah,  A.  J.  Berger,  D.  Jakubowski,  J.  B.  Fishkin,  J.  Butler,  and  B.  J.  Tromberg, 
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Reportable  Outcomes:  Additional  Funding 


We  have  received  funding  based  upon  this  work  from  the  California  Breast  Cancer  Research 
Program,  a  3  year,  $500,000  collaborative  effort  with  Dr.  Laura  Esserman  at  University  of 
California  at  San  Francisco. 

The  Chao  Comprehensive  Cancer  Center  has  also  received  a  gift  from  Avon  for  being  a  cancer 
center  for  excellence.  From  this  gift,  $250,0(X)  was  appropriated  to  the  optical  breast  cancer 
detection  program  because  of  our  contributions  to  the  field. 


Conclusions 

In  the  past  year,  we  have  made  some  significant  contributions  to  the  field  of  biomedical  optics 
that  may  ultimately  play  a  role  in  the  prevention,  detection,  and  treatment  of  breast  cancer.  Our 
work  provides  information  as  to  the  contrast  available  in  near-infrared  spectroscopy  of  the  breast. 
Our  results  on  normal  volunteers  demonstrate  adequate  physiologic  sensitivity  to  long-term 
hormone-induced  changes  in  the  breast.  We  have  just  begun  to  use  this  data  to  interpret  cancer 
measurements. 
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Near-infrared  light  offers  unique  possibilities  for  revealing  physiological  information 
from  within  biological  tissues  without  injections,  surgery,  dangerous  radiation,  or  great 
expense.  Although  it  has  been  known  for  some  time  that  non-ionizing  near-infrared  light 
penetrates  a  few  centimeters  into  biological  tissues,  only  recently  have  scientists  been 
able  to  take  advantage  of  this  by  taking  into  account  tissue  scattering.  Near-infrared  light 
provides  a  way  to  quantify  the  concentration  and  the  oxygenation  state  of  hemoglobin 
and  thus  characterize  breast  lesions  without  surgery.  This  information  will  provide 
clinicians  with  entirely  new  criteria  for  deciding  if  lesions  are  benign  or  malignant,  which 
may  help  reduce  some  of  the  ambiguities  and  false-negatives  currently  encountered  in 
conventional  breast  examinations. 

Optically  monitoring  physiological  changes  in  breast  tissue  also  provides  a  new  way  to 
look  at  both  disease  progression  and  prevention.  Using  near-infrared  light  from  an 
instrument  we  have  developed,  we  have  observed  changes  in  hemodynamics  and 
structure  that  are  characteristic  of  tumors.  In  addition,  we  have  also  observed  changes  in 
breast  tissue  physiology  resulting  from  hormone  replacement  therapy,  menopausal  status, 
and  estrous  cycle  stage.  The  technique  may  also  be  extended  to  monitoring  the  effects  of 
drugs  that  affect  breast  tissue  (such  as  estrogen  or  tamoxifen)  in  order  to  enhance  efficacy 
and  minimize  side  effects.  Because  our  technique  is  safe,  non-invasive  and  comparable 
in  cost  to  ultrasound,  patients  can  be  studied  as  often  as  necessary  without  risk. 
Currently,  we  are  performing  clinical  measurements  on  women  with  breast  tumors  and  on 
normal  volunteers  receiving  breast  cancer  prevention  drugs  such  as  tamoxifen  and 
raloxifene. 


The  U.  S.  Army  Medical  Research  and  Material  Command  under  DAMD 17-98- 1-8 186 
helped  support  this  work. 
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Near-infrared  (NIR)  optical  spectroscopy  demonstrates  unique  possibilities  as  a  non-invasive 
tool  for  monitoring  biological  tissue  physiology.  We  have  developed  a  bedside-capable 
instrument  that  uses  safe  levels  of  non-ionizing  NIR  light.  The  instrument  measures  the 
absorption  and  the  scattering  properties  of  tissues  by  applying  a  diffusive  model  of  light  transport 
to  photons  that  have  traveled  through  a  large  volume  of  tissue.  Hemoglobin,  water,  and  lipids 
may  be  quantified  by  their  absorption  properties,  and  scattering  information  is  sensitive  to  items 
such  as  cellular  density  and  collagen. 


We  have  used  a  prototype  instrument  uncovered  physiological  changes  in  breast  that  are 
consistent  with  known  changes  in  breast  parenchyma  as  a  result  of  age  and  menopausal  status  in 
a  sample  of  21  healthy  female  volunteers  (ages  18  to  60).  Combining  the  scattering  and  the 
absorption  information  together  reveals  functional  information  in  the  breast.  Such  information 
could  be  used  to  monitor  the  effectiveness  of  treatments  that  affect  the  composition  of  the  breast, 
such  as  hormone  replacement  therapy  and  neo-adjuvant  therapy.  All  of  this  was  done  with  an 
instrument  that  in  its  prototype  form  is  comparable  in  cost  to  commercialized  ultrasound 
technology. 
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Near-infrared  (NIR)  optical  spectroscopy  demonstrates  unique  possibilities  as  a  non-invasive 
monitor  of  tissue  physiology.  We  have  developed  a  bedside-capable  instrument  that  uses  safe 
levels  of  non-ionizing  NIR  light.  Our  instrument  measures  both  absorption  and  scattering 
properties  of  tissues  by  applying  a  diffusive  model  of  light  transport  to  photons  that  have 
traveled  through  a  large  volume  of  tissue.  Oxygenated  and  deoxygenated  hemoglobins,  water, 
and  lipids  were  quantified  by  their  absorption  signature.  Additional  information  was  obtained 
from  tissue  scattering,  which  is  sensitive  to  items  such  as  cellular  density  and  collagen  content. 


This  prototype  instrument  uncovered  physiological  changes  in  breast  that  are  consistent  with 
known  changes  in  breast  parenchyma  as  a  result  of  age  and  menopausal  status  in  a  sample  of  21 
healthy  female  volunteers  (ages  18  to  60).  Combining  scattering  and  absorption  information 
together  reveals  age-dependent  functional  information  in  the  breast.  This  physiological 
information  is  not  obtainable  with  any  other  method  to  date.  Such  information  could  be  used  to 
monitor  the  effectiveness  of  treatments  that  aflfect  breast  composition,  such  as  hormone 
replacement  therapy  and  neo -adjuvant  therapy,  as  well  as  evaluate  individual  breast  cancer  risk. 
The  prototype  instrument  in  its  present  form  is  comparable  in  cost  to  commercialized  ultrasound 
technology. 
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Abstract 

Frequency-domain  photon  migration  (FDPM)  is  a  non- 
invasive  optical  technique  that  utiiizes  intensity-modu¬ 
lated,  near- infrared  (NIR)  light  to  quantitativeiy  mea¬ 
sure  opticai  properties  in  thick  tissues.  Opticai 
properties  (absorption,  and  scattering,  ^5',  para¬ 
meters)  derived  from  FDPM  measurements  can  be  used 
to  construct  iow-resoiution  (0.5  to  1  cm)  functionai 
images  of  tissue  hemoglobin  (total,  oxy-,  and  deoxy- 
forms),  oxygen  saturation,  biood  voiume  fraction,  water 
content,  fat  content  and  ceiiuiar  structure.  Uniike 
conventionai  NIR  transillumination,  FDPM  enables 
quantitative  analysis  of  tissue  absorption  and  scattering 
parameters  in  a  singie  non -invasive  measurement.  The 
unique  functional  information  provided  by  FDPM  makes 
it  weil-suited  to  characterizing  tumors  in  thick  tissues. 
In  order  to  test  the  sensitivity  of  FDPM  for  cancer 
diagnosis,  we  have  initiated  clinicai  studies  to  quantita¬ 
tiveiy  determine  normai  and  maiignant  breast  tissue 
opticai  and  physioiogical  properties  in  human  subjects. 
Measurements  are  performed  using  a  non -invasive, 
muiti -wavelength,  diode-laser  FDPM  device  optimized 
for  clinical  studies.  Results  show  that  ductal  carcino¬ 
mas  (invasive  and  in  situ)  and  benign  fibroadenomas 
exhibit  1.25  to  3-fold  higher  absorption  than  normai 
breast  tissue.  Within  this  group,  absorption  is  greatest 
for  measurements  obtained  from  sites  of  invasive 
cancer.  Opticai  scattering  is  approximately  20%  greater 
in  pre-menopausai  versus  post -menopausal  subjects 
due  to  differences  in  gland/cell  proiiferation  and 
coliagen/fat  content.  Spatiai  variations  in  tissue  scatter¬ 
ing  reveal  the  loss  of  differentiation  associated  with 
breast  disease  progression.  Overall,  the  metabolic 
demands  of  hormonai  stimulation  and  tumor  growth 
are  detectabie  using  photon  migration  techniques. 
Measurements  provide  quantitative  opticai  property 
vaiues  that  refiect  changes  in  tissue  perfusion,  oxygen 
consumption,  and  ceil/matrix  deveiopment.  Neoplasia 
(2000)  2,  26-40. 

Keywords:  tissue  optical  properties,  absorption,  scattering,  diffuse  optical  imaging,  near- 
infrared  spectroscopy,  hemoglobin,  tumor  vasculature,  extracellular  matrix. 


Introduction 

Pathologists  routinely  examine  thin  sections  of  surgically 
removed  tissue  in  order  to  diagnose  cancer.  If  similar 
information  could  be  provided  by  techniques  that  allow  us 
to  see  internal  structure,  physicians  would  be  able  to  find 
and  characterize  tumors  non-invasively.  Historically,  this 
reasoning  lead  to  the  development  of  “diaphanography”,  a 
near-infrared  (NIR)  light  imaging  method  introduced  nearly 
70  years  ago  to  locate  and  identify  breast  cancer  [1].  NIR 
imaging  is  a  simple,  low-cost,  risk-free  procedure  that, 
unlike  X-ray  mammography,  does  not  employ  ionizing 
radiation.  Unfortunately,  the  initial  promise  of  NIR  diapha¬ 
nography  was  never  realized,  primarily  due  to  intense  light 
scattering. 

Distinguishing  between  absorption  and  scattering  is 
generally  recognized  as  one  of  the  principal  limitations  to 
detecting  cancer  with  light.  This  problem  has  challenged 
scientists  who  study  how  light  interacts  with  turbid  materials 
for  many  years.  Recent  developments  in  theory  and 
instrumentation,  collectively  referred  to  as  “photon  migra¬ 
tion”  (PM),  now  enable  us  to  probe  and  analyze  multiply- 
scattered  light  signals  on  an  ultra-fast  time  scale  (e.g., 
picoseconds).  These  technological  advances  provide  a 
simple  framework  for  measuring  the  exact  magnitude  of 
absorption  and  scattering  (i.e.,  optical  properties)  in  turbid 
materials.  Optical  properties  can,  in  turn,  be  used  to  locate 
and  identify  physiological  changes  characteristic  of  tumors, 
hormonai  stimulation,  pharmacologic  agents,  and  age- 
related  tissue  remodeling. 

Research  by  our  group  and  many  others  in  the  Biomedicai 
Optics  community  has  helped  rejuvenate  interest  in  opticai 
methods  for  detecting  and  characterizing  solid  tumors  [2]. 
The  use  of  ultra-fast  light  provides  new  information  with  the 
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same  advantages  as  conventional  NIR  diaphanography 
(e.g.,  low  cost,  zero  risk,  compact  devices).  However,  in 
order  to  demonstrate  whether  normal  and  malignant 
structures  can  actually  be  resolved  using  light,  reliable 
measurements  of  both  normal  and  tumor  tissue  optical 
properties  must  be  performed  in  human  subjects. 

Few  studies  have  been  completed  that  accurately 
measure  human  tissue  optical  properties  in  vivo.  This  can 
be  a  particularly  difficult  hurdle  to  overcome  because  photon 
migration,  as  a  fundamentally  new  method,  is  not  widely 
accessible.  PM  technology  development  requires  interdisci¬ 
plinary  expertise  in  diverse  fields  while  measurements 
require  strong  clinical  collaborations.  Clinical  studies  are 
now  underway  here  (about  70  human  subjects  have  been 
enrolled  in  breast,  brain,  and  cervical  protocols  to  date)  and 
in  other  laboratories  that  will  provide  critical  information  for 
determining  whether  non -invasive  photon  migration  mea¬ 
surements  are  useful  for  characterizing,  imaging,  and 
detecting  solid  tumors.  Early  results  show  that  unique  tumor 
optical  signatures  are,  in  fact,  measurable.  Organs  of 
greatest  initial  promise  appear  to  be  breast  and  brain. 
Changes  in  blood  flow,  oxygen  consumption,  and  tissue 
architecture  (cellular  and  matrix)  can  be  measured  in  these 


structures  quantitatively  using  NIR  photon  migration.  Most 
importantly,  these  signals  appear  to  reflect  fundamental 
physiological  processes  associated  with  malignancy,  dis¬ 
ease  progression,  and  response  to  therapy. 

In  this  work,  we  review  basic  principles  of  photon 
migration  [3]  in  thick  tissues  and  highlight  how  light-tissue 
interaction  mechanisms  provide  unique  opportunities  for 
solid  tumor  diagnostics.  Details  of  instrumentation  and 
theory  development  in  photon  migration  spectroscopy  are 
presented  along  with  results  of  clinical  studies  in  human 
breast  tumor  patients.  Because  many  attributes  related  to 
the  origin,  progression,  and  therapeutic  response  of  breast 
cancers  are  common  to  other  types  of  tumors,  these  data 
should  provide  insight  into  how  continued  optical  technology 
development  will  impact  cancer  diagnosis  and  screening. 


Light  Propagation  in  Thick  Tissues 

Spectrai  Dependence 

Light  propagation  in  biological  tissues  is  a  complex 
function  of  scattering  and  absorption,  which,  in  turn,  are 
dependent  on  cellular  structure  and  molecular  composition. 


Wavelength  (nm) 


Figure  1.  Upper  panels:  Spectral  dependence  of  tissue  chromophore  absorption.  A  tissue  optical  "window"  exists  in  the  red/ NIR  between  600  and  1300  nm  (from 
J.L.  Boulnois,  Lasers  in  Medical  Science,  1986: 1: 47-66).  Lower  panels:  Spectral  dependence  of  light  transmission  through  rabbit  ear  with  implanted  tumor  model. 
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Figure  1  illustrates  that  within  the  visible  (300  to  700  nm) 
water  window,  hemoglobin  and  melanin  are  the  dominant 
endogenous  absorbers  [4], 

Although  melanin  absorption  coefficients  are,  on  average, 
greater  than  those  of  hemoglobin,  the  overall  influence  of 
melanin  on  light  propagation  in  thick  tissues  is  generally  less 
substantial.  This  is  due  to  the  fact  that  melanin  is  typically 
confined  to  a  thin,  superficial  skin  layer.  In  contrast, 
hemoglobin  is  widely  distributed  throughout  most  tissues  at 
a  volume  fraction  that  ranges  from  1%  to  5%.  Exceptions  to 
this  include  avascular  structures  and  tissues  with  abnormally 
high  melanin  content,  such  as  melanotic  melanomas. 

Scattering  originates  from  inhomogeneities  in  tissue 
structure,  which,  in  turn,  are  determined  by  refractive 
index  discontinuities  occurring  both  between  and  within 
cells.  In  the  spectral  region  of  600  to  700  nm,  hemoglobin 
absorptivity  is  attenuated  by  approximately  20  dB  (nearly 
100-fold)  and  scattering  eclipses  absorption.  This  condi¬ 
tion  persists  until  water  absorption  resumes  as  the  primary 
attenuation  mechanism  at  approximately  1.3  /rm.  Thus, 
the  spectral  region  between  600  and  1300  nm  is 
considered  to  be  the  tissue  optical  “window”  since  both 
absorption  and  scattering  losses  are  minimal  throughout 
this  interval.  This  concept  is  illustrated  in  the  series  of 
three  photographs  shown  at  the  bottom  of  Figure  1.  Here, 
a  rabbit  ear  implanted  with  a  tumor  model  is  transillumi- 
nated  using  blue  (490  nm),  green  (530  nm)  and  red 
(630  nm)  light.  In  the  blue  and  green  photographs, 
excellent  contrast  is  observed  between  light-absorbing 
blood  vessels  and  the  surrounding  ear  cartilage.  This  is 
due  to  strong  hemoglobin  absorption  in  these  spectral 
regions  and  relatively  high  attenuation  of  light  from 
scattering.  Because  the  tumor  has  substantial  blood 
vessel/hemoglobin  content,  it  appears  as  a  large  light¬ 
absorbing  shadow.  In  the  red  light  case,  however,  the 
image  appears  quite  diffuse  and  the  small  vessel 
structures  have  disappeared.  This  is  a  consequence  of 
the  substantially  reduced  hemoglobin  absorption  at  630 
nm,  making  these  vessels  appear  to  be  “transparent”.  In 
addition,  red  light  is  scattered  less  than  blue/green.  As  a 
result,  many  of  the  scattered  photons  contribute  to  the 
image  as  “multiply  scattered  light,”  further  degrading 
spatial  resolution  and  contrast  while  improving  transpar¬ 
ency.  The  net  result  is  a  “cloudy”  tumor  that  appears  to 
be  different  than  its  actual  physical  dimension. 

The  wavelength -dependent  intensity  of  light  penetrance 
into  biological  tissues  can  be  used  to  characterize  structures 
of  different  composition.  Table  1  provides  several  practical 


examples  of  wavelength -dependent  “optical  penetration 
depths”,  6,  a  value  that  corresponds  to  the  distance  at  which 
the  optical  fluence  rate  (i.e.,  mW/cm^)  is  reduced  to  1/ 
e=0.37  of  the  initial  value  [5].  These  data  highlight  the 
dynamic  interplay  between  absorption  and  scattering.  For 
example,  since  human  retinoblastomas  contain  no  melanin, 
low  levels  of  hemoglobin,  and  are  only  minimally  scattering, 
penetration  depths  are  relatively  large;  ranging  from  2.9  mm 
at  600  nm  to  5.1  mm  at  1.06  /um.  Light  penetrance  is 
commensurately  reduced  in  brain  and  hand  tissue  due  to 
their  high-scattering/low-absorption  and  moderate -scatter¬ 
ing/moderate -absorption  properties,  respectively.  In  con¬ 
trast,  the  high  melanin  content  of  melanotic  melanomas 
strongly  attenuates  light  throughout  the  visible  region  and 
maximum  light  penetrance  is  only  1 .4  mm  at  1 .064  fim. 

These  examples  indicate  that,  in  general,  6  for  non- 
melanotic  tissues  follows  hemoglobin  absorption.  Thus,  from 
480  to  590  nm,  6  ranges  from  0.5  to  1.5  mm.  A  substantial 
penetration  depth  increase  from  about  1  to  3  mm  occurs 
between  600  and  650  nm,  and  a  more  gradual  increase  is 
observed  from  650  to  750  nm  (5=2  to  4  mm).  Between  750 
and  900  nm,  values  tend  to  plateau,  and  maximum  light 
penetrance  of  3  to  6  mm  generally  occurs  in  the  1 000  to  1 1 00 
nm  region. 

Multiple  Scattering 

On  a  cellular  level,  there  are  clear,  visible  differences 
between  normal  and  malignant  tissues.  These  changes  may 
include,  among  others,  multinucleation,  structural  distortion, 
hyperpigmentation,  and  membrane  and  mitochondrial  varia¬ 
tions.  On  a  bulk  tissue  level,  there  are  distinct  variations  in 
tumor  vasculature,  extracellular  matrix  and  biochemical 
milieu.  Pathologists  utilize  visible,  microscopic  anomalies  in 
cell/tissue  architecture  to  formulate  diagnoses.  Thus,  varia¬ 
tions  in  bulk  tissue  optical  penetration  depths  are  expected, 
given  the  large  body  of  supporting  histologic  evidence  [6]. 
These  fundamental  optical  property  differences  can,  in 
principle,  provide  novel  diagnostic  criteria  for  thick -tissue 
photon  migration  methods. 

However,  our  ability  to  employ  thick-tissue  optical 
methods  for  cancer  detection  is  limited  by  a  form  of  light¬ 
scattering  behavior  known  as  “multiple  scattering”.  Light 
launched  into  tissue  is  efficiently  scattered  out  of  an  incident 
collimated  beam  into  an  almost  isotropic  distribution  within  a 
few  millimeters  from  the  source.  This  “multiple  scattering” 
behavior  degrades  the  information  content  of  light  that 
propagates  through  thick  tissues.  In  practical  terms,  optical 
imaging  methods  that  cannot  distinguish  between  scattering 


Table  1.  Optical  Penetration  Depth  (mm)  as  a  Function  of  Wavelength  (From  Ref.  [5] ). 


Tissue 

Wavelength  (nm) 

600 

650 

700 

750 

800 

850 

900 

1064 

Human  retinoblastoma  {in  vitro) 

2.9 

3.8 

4.0 

4.0 

4.1 

4.2 

4.3 

5.1 

Porcine  brain  {in  vitro) 

1.8 

2.4 

2.9 

3,0 

3.3 

3.5 

3.7 

4.0 

Human  hand  {in  vivo) 

1.4 

2.0 

2.6 

2.7 

3.0 

3.0 

3.0 

_ 

Melanotic  melanoma  in  athymic  mice 

- 

0.28 

0.34 

0.41 

0.50 

0.56 

0.64 

1.4 
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and  absorption  are  incapable  of  providing  a  complete 
diagnostic  picture.  This  concept  is  illustrated  schematically 
In  Figure  2. 

In  homogeneous,  non -scattering  media,  light  propagation 

is  simply  a  function  of  the  absorption  coefficient,  and 

distance,  z: 

/=/oe-"^^  (1) 

where  i_ia=1/lat,  and  Ub  is  the  absorption  length,  or 
average  distance  between  absorption  events.  The  absorp¬ 
tion  coefficient  can  be  expressed  in  terms  of  molecular 
concentration  via  the  Lambert-Beer  law  (i.e.,  /4=log  II 
lo=EbC,  where  A,  e,  b,  are  C  are  absorbance  units,  molar 
extinction  coefficient,  pathlength,  and  molar  concentration, 
respectively): 

fi-a  =  2.3£C  (2) 

Thus,  measurements  of  light  penetrance  in  tissues  can,  in 
principle,  provide  quantitative  information  regarding  the 
concentration  of  physiologically  relevant  absorbers. 

However,  in  non-transparent,  turbid  media,  both  scatter¬ 
ing  and  absorption  contribute  to  the  distance -dependent 
light  attenuation.  The  total  attenuation  coefficient,  jUt,  is  the 
sum  of  absorption  and  scattering  parameters  (/(t=/<a+/-'s)- 
Here,  /j-s  is  the  scattering  parameter,  or  the  reciprocal  of  the 
scattering  length  (/sc).  In  the  special  case  of  “multiple 
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scattering”,  the  incident  light  intensity  diminishes  according 
to: 

I  =  (3) 

where  the  effective  attenuation  coefficient,  ^en='^l 
is  the  reciprocal  of  the  mean  penetra¬ 
tion  depth  (6).  In  order  to  account  for  multiple  scattering,  the 
reduced  scattering  parameter,  /u.s'=Ms(1  -0)  is  defined  in 
terms  of  //g  and  the  angular  dependence  of  scattering,  g, 
where  g={cos  theta),  (i.e,,  the  average  cosine  of  the 
scattering  angle).  In  most  tissues,  scattering  is  highly 
forward -directed  and  g  values  range  from  0.7  to  0.9  [6].  In 
practical  terms,  large  g  values  delay  the  distance  photons 
must  travel  before  optical  energy  is  isotropically  distributed. 
Generally,  this  “transport  mean  free  path”  (tmfp)  is  on  the 
order  of  1  mm  in  most  tissues  and  can  be  defined  as  the 
reciprocal  of  the  transport  scattering  coefficient,  ^tr=Ma+Ms'. 

The  Diffusion  Approximation 

In  order  to  probe  large  volumes  of  tissues  with  acceptable 
signal -to -noise  levels,  measurements  are  generally  con¬ 
ducted  in  spectral  regions  of  greatest  transparency,  i.e.,  the 
tissue  optical  window  ( 600  to  1 300  nml  see  Figure  1 ) .  Under 
these  conditions,  scattering  dominates  absorption  and  light 
propagation  can  be  described  by  a  diffusion  approximation 
[12,21],  Since  scatterers  and  absorbers  are  distributed  in 


Time  independent: 

Fluence  rate  decreases  exponentially 


Figure  2.  Red/NIR  light  launched  Into  tissue  propagates  with  directional  changes  (scattering)  until  captured  by  an  absorber  (molecule).  Multiple  scattering 
increases  the  mean  photon  path  length.  A  short,  e.g.,  picosecond,  light  pulse  launched  into  tissue  spreads  out  in  proportion  to  the  number  of  available  photon  paths. 
Absorption  is  a  loss  mechanism  that  reduces  the  number  of  available  photon  paths  and  hence  limits  the  temporal  dispersion  of  the  pulse.  Under  low  absorption 
conditions,  a  light  pulse  will  experience  maximal  dispersion.  The  photon  density,  or  fluence  rate  (ip)  decays  exponentially  with  distance  from  the  source.  The  1/e 
exponential  decay  constant,  or  optical  penetration  depth  (i)  is  a  function  of  the  absorption,  p  a,  and  scattering,  p  J,  properties  of  the  medium. 
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tissues  on  dramatically  different  spatial  scales  (i.e., 
/Us'z.Gt//a).  the  transport  scattering  length  can  be  approxi¬ 
mated  by  1  Under  these  conditions,  an  additional 

descriptive  domain,  i.e.,  time  or  frequency,  is  required  to 
resolve  light  scattering  from  absorption. 

As  shown  in  Figure  2,  launching  light  into  tissue 
results  in  multiple  scattering  events  that  can  dramatically 
increase  the  mean  path  length  beyond  the  linear  distance 
from  source  to  detector.  If  the  light  source  continuously 
illuminates  the  sample  (i.e.,  it  is  time -independent),  the 
fluence  rate  (in  mW/cm^)  or  number  density  of  photons 
measured  in  the  sample  from  all  directions  decays 
exponentially  with  distance  from  the  source.  The  decay 
rate  is  a  function  of  the  penetration  depth,  6,  which  in 
turn  can  be  given  by  the  absorption  and  scattering 
properties.  When  a  time -dependent  light  source  is  used 
such  as  a  picosecond  (psec)  laser  pulse  or  a  sinusoid¬ 
ally  intensity -modulated  laser,  the  large  number  of  photon 
paths  available  due  to  multiple  light  scattering  processes 
causes  broadening  of  the  source  temporal  profile.  Thus, 
measurements  of  the  time -dependent  propagation  char¬ 


acteristics  of  short  light  pulses  or  intensity- modulated 
waves  in  tissues  can  be  used  to  separate  the  contribu¬ 
tions  of  absorption  from  scattering.  When  combined  with 
an  appropriate  mathematical  framework,  this  allows 
precise  calculation  of  tissue  optical  properties,  namely 
tissue  absorption,  ^a,  and  reduced  scattering,  /ig', 
properties. 

Temporal  Resolution 

Over  the  past  few  years,  remarkable  gains  in  our 
understanding  of  tissue  optical  properties  have  been  realized 
by  interrogating  tissues  with  the  unique  temporal  (i.e.,  short 
pulse)  properties  of  lasers  [7-11].  Pulse  propagation 
methods  provide  information  about  the  distribution  of 
scatterers  and  absorbers  in  a  single  measurement  [12,13]. 
These  optical  properties  may  be  used  in  a  variety  of 
therapeutic  and  diagnostic  techniques,  including:  imaging 
tissue  structure  [14-16]  and  cortical  activity  [17,18]; 
monitoring  physiology  [19-23]  and  drug  concentration 
[24];  dysplastic  transformation  [25],  and  predicting  optical 
dosimetry  for  laser-based  procedures  [26]. 


Source  Detector 
I  25  mm  I 

I  i 


tj  =  200  ps  t2  =  600  ps  t3  =  1000  ps  =  1400  ps 


Figure  3.  Numerical  simulation  showing  time -dependent  propagation  of  ultrashort  light  pulse  in  tissue  with  optical  properties  similar  to  human  breast 

cm^  \  p  s'  =10cm~  ^ ).  Two-dimensional  x-zviews  (40x25  mm^)  of  detected  photon  density  are  shown.  Time -gating  of  the  detector  provides  control  over  depth 

and  volume  of  tissue  probed. 
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The  conceptual  basis  for  the  time -domain  approach 
generally  involves  solutions  to  the  radiative  transfer 
equation  [27,28]  using  Monte  Carlo  simulation  [29-31] 
and  diffusion  theory  approximations  [7,32].  Diffusion - 
based  methods  provide  relatively  straightforward  analytical 
expressions  that  describe  the  shape  of  a  diffusely  reflected 
or  transmitted  pulse  in  terms  of  the  optical  properties  of  the 
medium  [10].  Thus,  the  observed  temporal  broadening  of 
ultra  short  pulses  can  be  mathematically  related  to  the 
large  number  of  optical  paths  available  in  multiple - 
scattering  media.  Since  the  introduction  of  losses  (absor¬ 
bers)  reduces  the  average  path  length,  absorber-depen- 
dent  changes  in  pulse  propagation  time  can  be  used  to 
calculate  absorption  coefficients  [10].  This  concept  is 
illustrated  in  Figure  3.  Here,  a  picosecond  light  pulse  is 
launched  into  a  tissue -like  medium  with  optical  coefficients 
similar  to  human  breast  tissue  (/ta=0.1  cm“\  ^is'=10 
cm“^).  The  transport  characteristics  of  the  source  are 
simulated  numerically  using  a  Monte  Carlo  technique  [33]. 
The  figure  shows  how  the  detected  pulse  spreads  out  in 
time  and  space,  due  to  the  large  number  of  photon  paths 
available.  By  time -gating  the  detector,  one  can  adjust  the 
“viewing  volume”  in  the  tissue.  For  example,  at  early  times 
after  the  pulse  is  launched,  shorter  photon  paths  are 
detected.  This  is  demonstrated  by  the  200  psec  time  point 
that  samples  tissue  volumes  above  the  1  cm  depth  line. 
When  the  200  psec  observation  “gate”  is  shifted  to  f=600 
psec,  the  mean  interrogation  depth  and  volume  are 
substantially  greater.  This  is  due  to  the  fact  that  at  longer 
times,  photons  have  traveled  tens  of  centimeters,  subs¬ 
tantially  further  than  the  2.5  cm  linear  distance  between 
source  and  detector.  By  moving  the  observation  window 
through  the  time  course  of  the  detected  signal,  the  spatial 
distribution  of  absorbers  and  scatterers  in  the  tissue  can  be 
reconstructed,  and  sensitivity  to  deep  tissue  structures, 
such  as  tumors,  can  be  optimized. 

Frequency-domain  optical  methods  can  be  adapted  to 
diffusion  theory  models  in  a  similar  manner.  Fishkin  and 
Gratton  [34]  first  suggested  that  amplitude -modulated  light 
propagates  through  homogeneous  multiple  -  scattering  media 
as  diffuse  waves  with  a  coherent  front.  These  photon  density 
waves  can  be  characterized  by  a  phase  velocity  (Vp)  and 
modulation  wavelength  (Am)  that  are  primarily  functions  of 
media  optical  properties.  Measurements  are  performed  by 
launching  intensity -modulated  light  into  tissue  and  recording 
the  phase  delay  {phi)  and  demodulation  amplitude  (m)(with 
respect  to  the  source )  at  a  fixed  distance  from  the  launch  site. 
Optical  properties  (the  absorption  coefficient,  and  the 
reduced  scattering  coefficient,  ps)  are  calculated  from  the 
measured  frequency-  and  distance -dependent  behavior  by 
fitting  the  experimental  results  to  an  appropriate  theoretical 
framework.  Thus,  FDPM  can  rapidly  and  quantitatively 
assess  tissue  optical  properties  in  a  single  non -invasive 
measurement.  It  is  important  to  point  out  that  time  and 
frequency  domain  methods  are  equivalent  and  related  to  each 
other  through  the  Fourier  transform.  We  employ  the 
frequency-domain  approach  due  to  instrumentation  advan¬ 
tages  that  include  reduced  cost,  complexity,  and  simplicity  of 


mathematical  techniques  used  in  frequency-domain  signal 
analysis  [35]. 


Instrumentation 

In  frequency -domain  photon  migration  (FDPM),  the  inten¬ 
sity  of  light  incident  on  an  optically  turbid  sample  is 
modulated  at  high  frequencies  (e.g.,  hundreds  of  mega¬ 
hertz),  and  the  diffusely  reflected,  transmitted,  or  re -emitted 
(e.g.,  fluorescent)  signal  is  measured  with  a  phase - 
sensitive  detector.  Intensity -modulated  light  propagates 
through  multiple -scattering  media  with  a  coherent  front, 
forming  photon  density  waves  (PDWs).  PDW  dispersion  is 
highly  dependent  on  the  optical  properties  of  the  medium. 
Exact  absorption  (pa)  and  reduced  scattering  (ps)  para¬ 
meters  are  calculated  by  comparing  the  measured  fre¬ 
quency-  or  distance -dependent  PDW  phase  and  amplitude 
behavior  to  analytically  derived  non-linear  model  functions. 
Model  functions  are,  in  turn,  derived  from  solutions  to  the 
photon  diffusion  equation  which  yield  expressions  for  phase 
(phi)  and  amplitude  (A)  as  a  function  of  modulation 
frequency  (w)  and  tissue  optical  properties  under  various 
boundary  conditions  [36].  The  values  phi,  A,  and  u  are 
defined  in  the  usual  manner  where  ph/=the  measured  phase 
lag  between  the  source  (reference)  and  sample  response; 
/A=/4Csampie/7tCsouroe;  and  w=27rf  is  the  angular  modulation 
frequency. 

A  schematic  of  our  1-GHz,  portable  FDPM  device  is 
shown  in  Figure  4  [37].  FDPM  instrumentation  and 
theoretical  background  have  been  described  in  detail  else¬ 
where  (Fishkin  et  al.  [41];  Chance  et  al.  [35]).  Briefly, 
modulation  is  swept  from  300  kHz  to  1  GHz  in  less  than  1 
second,  providing  rapid,  multi -wavelength  (~670  to  960 
nm)  characterization  of  most  tissues  in  a  single  measure¬ 
ment.  The  instrument  is  compact  and  can  easily  be 
transported  to  operating  rooms  and  bedridden  patients.  By 
incorporating  several  NIR  laser  diodes,  a  number  of 
wavelength -dependent  physiological  parameters  are  deter¬ 
mined.  Source  modulation  frequencies  typically  range  from 
300  kHz  to  1  GHz,  with  FDPM  data  recorded  in  5  MHz 
increments. 

The  core  component  of  the  FDPM  apparatus  is  a  network 
analyzer  (Hewlett  Packard,  model  8753C),  which  is  used  to 
produce  modulation  swept  from  300  kHz  to  1  GHz  (20  dB  m 
RF  output).  RF  from  the  network  analyzer  is  serially 
superimposed  (via  the  AC  switch)  on  the  direct  current  of 
up  to  eight  different  diode  lasers  (e.g.,  SDL,  Inc.  models 
7421 , 5420,  5421 ,  and  6321  at  674,  81 1 , 849,  and  956  nm, 
respectively)  using  individual  bias-tees  (model  5575A, 
Picosecond  Pulse  Labs)  and  an  RF  switch  (model  8768K, 
Hewlett  Packard).  One -hundred -micrometer- diameter  gra¬ 
dient-index  fibers  are  used  to  couple  each  light  source  to  an 
8x8 optical  multiplexer  (model  GP700,  DiCon  Instruments). 
The  8x8  optical  multiplexer  allows  for  up  to  eight  different 
diode  laser  light  sources  and  eight  different  optical  fiber 
positions. 

Light  is  launched  onto  the  tissue  (or  test  object)  using  the 
above-mentioned  unique  wavelengths  and  one  source  fiber. 
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Figure  4.  One -gigahertz,  multi -frequency,  multi  -  wavelength  FDPM  Instrument.  See  text  for  complete  description. 


An  avalanche  photodiode  (APD,  Hamamatsu,  model 
C5658)  is  used  to  detect  the  diffuse  optical  signal  that 
propagates  through  the  biological  tissue.  Both  the  APD  and 
probe  end  of  the  source  optical  fiber  are  in  direct  contact  with 
the  patient  (i.e.,  a  “semi -infinite  medium”  measurement 
geometry).  The  optical  power  coupled  into  the  tissue 
averages  approximately  10  to  30  mW.  Measurement  time 
depends  on  the  precision  required,  the  number  of  sweeps 
performed,  and  RF/optical  switch  times.  For  human  subject 
studies,  approximately  0.1  second  is  used  to  sweep  over  the 
entire  1  GHz  band  of  modulation  frequencies.  However,  total 
elapsed  time  for  four  diodes  (typically  12  to  16  sweeps/ 
diode),  data  transfer,  display,  and  source  switching  is 
approximately  40  seconds.  Most  components,  including  the 
network  analyzer,  RF/optical  switches,  diode  power  supplies 
and  temperature  of  diode  mounts  are  controlled  by  computer 
using  virtual  instrument  software  (LabView,  National  Instru¬ 
ments). 

As  shown  schematically  in  Figure  5,  phase  and 
amplitude  data  (represented  by  phi  and  A,  respectively) 
obtained  from  frequency-dependent  FDPM  tissue  mea¬ 
surements  are  fit  to  theoretical  model  functions.  This  results 
in  calculation  of  the  optical  absorption  coefficient,  and 
the  reduced  scattering  coefficient,  ps\  at  a  given  A  and 
source-detector  separation  p.  Minimization  of  the 
surface  is  achieved  for  simultaneous,  error-weighted  fitting 
of  the  phi  and  A  versus  frequency  data  using  a  Marquardt- 


Levenberg  algorithm.  Typical  pa:  and  ps  uncertainties, 
determined  from  the  distribution  of  phase  and  amplitude 
fits,  range  from  0.5%  to  5%  of  the  mean  value.  When  optical 
properties  {ps  and  pa)  are  recovered  for  various  source 
wavelengths,  the  spectral  dependence  of  absorption  can  be 
combined  with  known  values  of  molecular  extinction 
coefficients  to  calculate  physiologically  relevant  parameters, 
such  as  oxygenated,  deoxygenated,  and  total  hemoglobin 
concentration;  oxygen  saturation;  drug  concentration;  blood 
volume  fraction;  and  water  concentration.  The  spectral 
dependence  of  scattering  is  used  to  provide  insight  into  the 
mean  density  and  size  of  scattering  structures. 


Application  to  Solid  Tumor  Characterization 

The  sensitivity  of  optical  properties  to  tissue  structure  and 
function  suggests  that  differences  between  normal  and 
diseased  tissues  may  provide  a  means  for  assessing  the 
probability  that  processes  such  as  malignant  transformation, 
inflammation,  or  infection  are  occurring.  Indeed,  previous 
studies  have  shown  unique  scattering  and/or  absorption 
signatures  associated  with  dysplastic  [25,38,39],  malignant 
[41]  and  benign  [45]  tissue  transformations. 

Characterization  of  solid  tumors  in  thick  tissues  poses 
unique  challenges.  These  structures  are  typically  found 
embedded  beneath  epithelial  layers  and  require  methods 
capable  of  probing  deep  tissue  structure  and  function. 
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Experimental  Response 
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Figure  5.  FDPM  measurement  overview:  intensity -modulated  light  is  launched  into  tissue,  resulting  in  the  propagation  of  diffuse  photon  density  waves.  The 
frequency -dependent  experimental  response  (phase,  phi,  and  amplitude,  A  i/s.  frequency,  lo)  is  compared  to  a  theoretical  model.  Model -data  fits  provide 
calculated  values  for  absorption,  /i  a,  and  reduced  scattering,  //  J,  coefficients  at  each  optical  wavelength.  Absorption  and  scattering  spectra  are  generated. 
Absorption  spectra  are  used  to  calculate  concentration  of  principal  tissue  chromophores. 


Because  of  multiple  light  scattering,  spatial  resolution  using 
“diffuse  methods”  is  intrinsically  rather  poor,  on  the  order 
of  1  to  5  mm  [40].  Nevertheless,  Imposing  temporal 
resolution  can  confine  the  photon  Interrogation  volume  to 
specific  regions  (see  Figure  3)  and  improve  sensitivity  to 
buried  tumors.  This  is  achieved  by  examining  the  temporal 
point  spread  function  of  a  short  light  pulse,  as  Illustrated  In 
Figure  3,  or  by  measuring  the  phase  and  amplitude  of 
photon  density  waves,  as  illustrated  In  Figure  5.  By 
measuring  phi  and  A  at  multiple  modulation  frequencies 
(w),  information  that  is  comparable  to  recording  the 
temporal  profile  of  the  diffusely  propagating  pulse  is 
obtained.  When  these  data  are  recorded  at  multiple 
wavelengths  (A),  FDPM  becomes  a  quantitative  spectro¬ 
scopic  technique. 

Structural  components  that  contribute  to  FDPM  signals 
in  thick  tissues  are  shown  in  Figure  6-4.  The  appearance 
and  progression  of  tumors  alter  the  fractional  contribution 
of  cells,  extracellular  matrix,  and  vasculature  to  the 
absorption  and  scattering  spectra.  In  addition,  spatial 
variations  In  optical  properties  may  be  evident  for  poorly 


differentiated  tumor  versus  well -differentiated  normal  tis¬ 
sues.  This  is  due  to  the  fact  that,  with  disease  progression, 
there  is  typically  an  Increase  In  tumor  cell  density,  nuclear 
volume  fraction,  and  stromal/inflammatory  cells.  Extracel¬ 
lular  matrix  degrades  with  increasing  tumor  invasion,  and 
neovasculature  Is  recruited,  forming  a  tortuous,  leaky 
network  complete  with  flow  anomolles.  These  alterations 
in  tissue  architecture  affect  both  the  absolute  value  and 
wavelength  dependence  of  the  reduced  scattering  para¬ 
meter.  From  a  molecular  perspective,  we  assume  that  the 
chromophores  contributing  to  at  a  given  NIR  wave¬ 
length  are  principally  oxy-  and  deoxy- hemoglobin,  water, 
and  fat.  The  total  absorption  measured  at  a  given 
wavelength  is  a  linear  sum  of  each.  Thus,  in  order  to 
achieve  optimal  sensitivity,  it  is  desirable  to  illuminate  the 
tissue  with  light  source  wavelengths  that  correspond  to 
high-contrast  absorption  features.  For  example,  in  Figure 
66,  fractional  spectral  contributions  are  shown  for  physio¬ 
logical  concentrations  of  principal  components  In  a  typical 
breast  tissue:  deoxy-  and  oxy-hemoglobin,  fat,  and  water. 
Sources  at  about  670,  810,  920,  and  960  nm  are  most 
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Figure  6.  fA)  Structural  and  functional  components  of  the  optical  signal  in  normal  and  transformed  tissue.  (B)  Absorption  spectra  of  tissue  with  composition  similar- 
to  human  post  -  menopausal  breast:  10  pM  deoxyhemoglobin  [Hb],  15  pM  oxyhemoglobin  [Hb02],  fat  (80%  volume  fraction ) ,  and  water  (15%  volume  fraction ) . 


sensitive  to  each.  In  order  to  calculate  tissue  concentra¬ 
tions,  Pa  measurements  at  a  minimum  of  four  different 
wavelengths  are  combined  with  known  values  of  molecular 
extinction  coefficients  [41  ].  Due  to  the  metabolic  demands 


of  rapidly  proliferating  tumor  cells  and  increased  tumor 
vasculature,  these  biochemical  constituents  can  vary 
substantially  in  relative  abundance  for  normal  versus 
diseased  tissue. 
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Breast  Measurements:  Results  and  Discussion 

The  relative  impact  of  these  physiological  features  on 
photon  migration  can  be  determined  by  examining  normal 
breast  tissue  optical  properties  (see  Ref.  [42]  for  a 
description  of  breast  tissue  optics).  Figure  7  shotws  results 
of  FDPM  measurements  obtained  from  two  pre- menopausal 
and  two  post- menopausal  subjects  (27, 29  and  63, 67  years 
old,  respectively).  Data  were  collected  from  the  right  upper 
outer  quadrants  of  each  subject  while  in  a  reclining  position. 
The  probe  was  applied  with  gentle  pressure  (i.e.,  no  external 
compression)  using  a  source-detector  separation  of  2.5 
cm. 

Pre-  and  post -menopausal  normal  breasts  exhibit  clear 
wavelength -dependent  differences  in  both  absorption  and 
scattering  parameters.  Figure  76  shows  that  scattering  is 
approximately  20%  greater  in  pre -menopausal  (//s'~0.8  to 
1.1  mm“^)  versus  post- menopausal  women  (^(s'~0.6  to 
0.7  mm“"^).  In  addition,  the  wavelength  dependence  of 
scattering  is  steeper  for  pre-  versus  post -menopausal 
subjects.  These  findings  are  consistent  with  known  changes 
in  breast  physiology.  Breast  tissue,  while  under  hormonal 
control,  has  higher  glandular/cellular  content  and  collagen  is 
required  in  the  extracellular  matrix  in  order  to  support  the 
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Figure  7.  FDPM  measurements  in  PRE  (9)  and  POST  (O)  menopausal 
normal  breast  (four  subjects).  (A)  Calculated  values  of  absorption 
coefficients,  /i  a,  vs.  wavelength.  (B)  Calculated  values  of  reduced  scattering 
coefficients,  ps,  vs.  wavelength. 


Structural  demands  of  cellular  proliferation.  After  meno¬ 
pause,  the  absence  of  hormonal  stimulation  results  in 
glandular  shrinkage,  collagen  remodeling  to  fat,  and  dimin¬ 
ished  vasculature.  Since  small  tissue  structures  such  as 
collagen  fibers  and  subcellular  organelles  are  likely  to  be  the 
primary  contributors  to  tissue  scattering,  it  is  not  surprising 
that  the  wavelength  dependence  of  scattering  is  relatively 
steep  for  pre -menopausal  subjects.  In  contrast,  following 
conversion  to  a  principally  large -particle  fatty  matrix,  a  much 
more  gradual  wavelength  dependence  of  scattering  is 
observed.  Thus,  the  diminished  contribution  of  the  glands 
and  collagenous  stroma  leads  to  a  clear  wavelength - 
dependent  scattering  reduction. 

Variations  in  absorption  can  be  explained  in  a  similar 
manner.  Absorption  contrast  due  to  enhanced  perfusion  is 
clearly  detectable  while  breast  tissue  is  under  hormonal 
control.  This  is  reflected  in  Figure  7A  values  of  absorption 
coefficients  that  are  consistently  higher  for  pre -menopausal 
subjects.  Following  menopause,  the  absence  of  substantial 
metabolic  demand  leads  to  reduced  perfusion  and  lower 
absorption  values.  This  is  well -correlated  with  the  low  blood, 
high  fat  content  of  post- menopausal  tissue.  Similar  changes 
in  optical  properties  have  been  observed  by  others  [43,44] . 

These  results  demonstrate  the  intrinsic  sensitivity  of 
FDPM  to  breast  tissue  physiological  states.  In  order  to 
address  the  origin  of  photon  migration  signals  in  solid 
tumors,  we  have  initiated  breast  tumor  studies  in  human 
subjects  [45] .  Our  goal  is  to  carefully  define  the  level  of 
accuracy  and  precision  required  from  our  measurement 
techniques  so  that  malignant  lesions  can  be  identified  and 
resolved  from  normal  structures  with  a  reasonable  level  of 
confidence.  The  practical  impact  of  increasing  diagnostic 
specificity  would  be  two -fold:  1 )  enhanced  early  detection  of 
relatively  small  (0.5to1  cm  diameter)  tumors  is  likely  to  lead 
to  mortality  improvements  (mammographic  performance  is 
highly  variable  in  pre -menopausal  subjects);  and  2)  more 
specific  characterization  of  benign  versus  malignant  lesions 
would  lead  to  reductions  in  the  large  number  of  unnecessary 
surgical  biopsies.  As  a  result,  there  is  considerable  room  for 
the  development  of  new,  non -invasive  optical  methods  for 
characterizing  breast  tissue. 

In  the  following  discussion,  we  present  results  from 
three  individuals.  Experiments  were  performed  under  the 
guidelines  of  UC  Irvine  I RB- approved  protocol  #95-563. 
Patient  1  was  a  67 -year -old  post -menopausal  subject 
with  a  single  palpable  mass  7.4  mm  beneath  the  skin 
surface  in  the  upper  outer  quadrant  of  the  left  breast. 
Histological  examination  following  surgical  biopsy  and 
ultrasound  revealed  a  1.8x0. 9  cm  ductal  carcinoma  in 
situ  (DCIS)  (malignant  tumor).  Patient  2  was  a  63-year- 
old  post -menopausal  subject  with  a  non-palpable  spicu- 
lated  left  breast  mass.  Histological  examination  following 
surgical  biopsy  revealed  a  0.6x0. 4x0. 7  cm  invasive  ductal 
carcinoma  (malignant  tumor).  Patient  3  was  a  22-year- 
old,  pre- menopausal  subject  patient  with  a  2x2  cm 
palpable  mass  3-7  mm  below  the  surface  in  the  upper 
medial  portion  of  the  right  breast.  Surgical  biopsy  revealed 
the  presence  of  a  benign  fibroadenoma. 
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Figure  8.  FDPM  measurements  of  phase  (A)  and  amplitude  (B)  vs.  source 
modulation  frequency  obtained  from  normal  and  tumor  sites  on  patient  1 
(ductal  carcinoma  in  situ,  DCIS ) .  Source-detector  distance =2.5  cm;  A = 674 
nm.  Solid  lines  are  best  simultaneous  fits  to  phase  and  amplitude  data  from 
tumor  (*)  and  normal  (O)  data  points,  respectively.  Ten  percent  of  200 
frequency -domain  data  points  is  shown  for  clarity. 


Measurements  are  performed  on  each  patient  by  gently 
placing  the  FDPM  probe  on  both  normal  and  tumor- 
containing  breast.  Data  are  acquired  using  a  hand-held 
scanning  probe  placed  in  nine  discrete  locations  covering 
a  2x2  cm®  grid  mapping  the  breast  surface.  The  probe 
incorporates  source-detector  (s-d)  pairs  ranging  from  1 
to  2.5  cm  in  separation.  The  2.5  cm  s-d  configuration  is 
placed  on  the  tissue  with  the  source  and  detector 
bracketing  the  tumor.  Tumor  location,  dimension,  and 
depth  are  monitored  immediately  prior  to  FDPM  measure¬ 
ments  by  ultrasound.  Photon  migration  data  are  acquired 
by  moving  the  probe  in  0.5  mm  increments  along  inferior- 
superior  and  medial -lateral  paths.  Repeat  measurements 
immediately  above  tumor  center  are  obtained  at  least 
three  times.  Both  normal  and  tumor-containing  breasts 
are  studied.  Each  tissue  location  is  characterized  in  terms 
of  wavelength -dependent  absorption  and  scattering  para¬ 
meters  which,  in  turn,  are  used  to  calculate  physiological 
properties. 

Sequential  scans  of  the  same  location  following  probe 
removal  and  replacement  revealed  no  significant  variations 
in  optical  properties.  Normal  tissue  measurements  are 


acquired  in  the  same  manner  from  a  symmetric  site  on  the 
opposite,  uninvolved  breast.  Phase  and  amplitude  data 
(represented  by  phi  and  A,  respectively)  obtained  from 
FDPM  tissue  measurements  are  compared  to  semi -infinite 
model  functions  to  extract  the  optical  absorption  coefficient, 
Pa,  and  the  reduced  scattering  coefficient,  ps,  at  a  given  A 
and  source-detector  separation.  Typical  pa  and  ps  un¬ 
certainties,  determined  from  the  phase  and  amplitude  fits, 
range  from  0.5%  to  5%  of  the  mean  value. 

Results  of  674-nm  FDPM  measurements  are  shown 
for  patient  1  in  Figure  8.  Raw  data  reveal  clear  differences 
in  both  phase  (Figure  8/4)  and  amplitude  (Figure  8B)  for 
normal  and  tumor  sites.  Solid  lines  are  results  of  model- 
fits  data  demonstrating  good  agreement  between  mea¬ 
surements  and  the  photon  migration  theoretical  model. 
Only  10%  of  the  actual  data  points  is  shown  in  each 
figure  for  clarity. 

Optical  properties  calculated  from  each  scan  position  are 
compared  to  mean  values  obtained  from  11  discrete 
locations  on  the  normal  side.  Figure  9/4  presents  the  ratio 
of  absorption  coefficients  (tumor/normal)  acquired  in 
medial -lateral  and  superior-inferior  scans  at  674  nm. 
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Figure  9.  (A)  Ratio  of  absorption  coefficients  ( tumor/ normal )  acquired  in 
medial -lateral  and  superior- inferior  scans  at  674  nm  for  patient  1  (ductal 
carcinoma  in  situ,  DCIS ) .  Normal  side  represents  average  of  11  mulfiple  site 
measurements.  Contrast  of  approximately  three-  fold  is  observed  for  tumor  vs. 
normal  sites  when  the  probe  is  placed  5  mm  lateral  of  the  estimated  tumor 
center.  (B)  Wavelength  dependence  of  optical  properties  acquired  from  the 
peak  contrast  location  (5  mm  lateral,  medial -lateral  axis)  and  symmetric 
location  on  normal  side. 
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Contrast  of  approximately  three-fold  is  observed  for  tumor 
versus  normal  sites  when  the  probe  is  placed  just  5  mm 
lateral  of  the  tumor  center.  A  clearly  discernable  tumor  shape 
is  apparent  in  this  horizontal  scan  direction.  Reduced  tumor/ 
normal  absorption  contrast  is  observed  in  the  superior- 


inferior  scan,  most  probably  due  to  differences  in  tumor 
orientation  along  this  axis. 

Wavelength -dependent  optical  properties  acquired  from 
the  peak  contrast  location  (5  mm  lateral,  medial -lateral 
axis)  are  displayed  in  Figure  9S.  Absolute  tumor 
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Figure  10.  (A)  Hemoglobin  concentration  (^iM)  in  deoxy-  [Hb],  oxy-  [HbOg],  and  totai  forms;  Xoxygen  saturation  ( 100x  [Hb02]/ IHb(tot)  ] ),  and  water 
content  (%)  for  normal  and  tumor  breast,  patienlt .  Values  calculated  from  wavelength  dependence  of  absorption  (Figure9B).  (B)  The  ratio  of  p  J  BSD  values  for 
normal /tumor  tissues  at  each  measurement  wavelength  for  patients  1  ( DCIS ),  2  ( invasive  ductal  carcinoma ) ,  and  3  (benign  fibroadenoma ) .  BSD  is  determined 
from  the  mean  and  standard  deviation  of  11  measurements  and  nine  discrete  locations  on  normal  and  tumor  sides.  A  BSD  ratio  of  1  implies  that  tumor  and  normal 
tissue  have  equivalent  structural  variation. 
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values  range  from  a  minimum  of  about  0.007  mm^^  at 
803  nm  to  a  maximum  of  0.017  mm'^  at  956  nm. 
Regardless  of  wavelength,  tumor  absorption  is  consis¬ 
tently  on  the  order  of  2-  to  3.5-fold  greater  than  normal 
tissue.  This  is  due  to  the  enhanced  tumor  versus  normal 
blood  supply  (perfusion)  and  higher  hemoglobin  content. 
Tumor  and  normal  tissue  absorption  spectra  also  have 
different  shapes.  The  slope  between  674  and  803  nm 
features  is  negative  for  tumor,  positive  for  normal  tissue. 
The  significance  of  this  is  related  to  the  increased 
contribution  of  de- oxyhemoglobin  to  light  absorption  in 
spectral  regions  that  are  shorter  in  wavelength  than  the 
~800  nm  isosbestic  point  (i.e.,  the  point  of  equivalent 
absorption  efficiency  for  both  oxy-  and  deoxy- hemoglobin 
forms). 

Physiologic  differences  between  tumor  and  normal  tissue 
are  further  amplified  by  calculations  of  water  and  blood 
content.  The  concentrations  of  water,  deoxy-,  oxy-,  and 
total  hemoglobin  in  DCIS  and  normal  breast  are  determined 
from  known  extinction  coefficients  and  Figure  96  absorption 
data.  Results,  shown  in  Figure  10/4,  clearly  reveal  elevated 
tumor  levels  for  each  physiological  parameter.  Total  tumor 
hemoglobin  ( [Hb(tot)  ]  =  [Hb]  +  [Hb02] )  is  35  ^M,  nearly 
twice  the  20  level  typical  of  normal  postmenopausal 
breast.  Deoxy- hemoglobin  ([Hb])  values  of  about  15  pM 
are  roughly  three-fold  greater  than  normal.  Together,  these 
values  can  be  used  to  determine  percent  oxygen  saturation, 
Y,  where  V%  =  100x[HbO2]/[Hb(tot)].  In  this  case, 
/(tumor) ~65%  and  /( normal )~80%.  Differences  are 
due  to  the  oxygen  extraction  demands  of  rapidly  proliferat¬ 
ing,  metabolically  active  tumor  cells.  The  reduced  /value  is 
consistent  with  the  measured  elevation  in  hemoglobin  at  the 
tumor  since  both  blood  and  oxygen  are  required  to  sustain 
tumor  growth.  Others  have  reported  similar  changes  in 
tumor  hemoglobin  concentration  and  oxygen  saturation 
using  non -invasive  photon  migration  techniques  in  the 
breast  [46,47]. 

Tissue  water  concentration  is  also  displayed  in  Figure 
10/4.  It  is  difficult  to  confirm  the  accuracy  of  these  values 
since  they  are  based  on  pure  water  extinction  coefficients 
(as  opposed  to  protein -bound  forms)  at  25°C  and  our 
calculations  do  not  take  into  account  the  contribution  of  fat  to 
the  956  nm  signal.  Nevertheless,  results  fall  within  the  1 1 .4% 
to  30.5%  range  given  in  the  literature  for  percentage  water  in 
human  fatty  adipose  tissue  [48] .  In  addition,  tumors  develop 
high  interstitial  fluid  pressure  due  to  lack  of  lymphatic 
drainage.  Thus,  increased  fluid  retention,  blood  content, 
and  cellularity  are  likely  to  be  the  cause  of  the  slightly 
elevated  tumor  water  percentage  observed  in  this  patient 
(20%  vs.  15%). 

Despite  these  clear  biochemical  differences,  it  is  not  yet 
known  whether  malignant  tumors  can  be  distinguished  from 
benign  lesions  using  this  technology.  In  fact,  it  is  our 
experience  that  both  malignant  and  benign  tumors  generally 
display  elevated  hemoglobin  and  reduced  /%  values. 
Although  these  observations  are  for  a  limited  number  of 
patients,  it  is  likely  that  additional  diagnostic  specificity  is 
required.  This  can  be  provided  by  examining  spatial 


variations  in  the  tissue -scattering  parameter,  fis-  As 
illustrated  schematically  in  Figure  8,  normal  well -differen¬ 
tiated  tissue  is  composed  of  distinct  structural/functional 
domains.  The  appearance  of  disease  perturbs  this  intrinsic 
organization,  resulting  in  a  gradual  loss  of  structural 
heterogeneity.  Ultimately,  the  growth  of  an  invasive  tumor 
will  lead  to  a  chaotic,  poorly  differentiated  environment. 

The  precise  relationship  between  disease  progression 
stage  and  photon  migration -scattering  signature  is  poorly 
understood.  In  order  to  explore  this  more  carefully,  we  have 
mapped  tumor  optical  properties  from  three  patients  with 
benign,  malignant  invasive,  and  malignant  in  situ  breast 
disease.  As  described  above,  data  are  acquired  from  at  least 
11  measurements  (nine  discrete  locations  in  5  mm  incre¬ 
ments)  using  a  2x2  cm^  scanning  pattern  on  the  breast 
surface.  Absorption  and  scattering  coefficients  are  com¬ 
puted  and  mean  ± standard  deviation  (SD)  values  calcu¬ 
lated  for  all  measurements  of  normal  and  tumor- 
containing  sites.  The  relative  standard  deviation  (RSD) 
of  optical  properties  obtained  from  tumor  and  normal  sites 
is  expressed  in  the  usual  manner,  RSD=100x(SD/ 
mean).  Figure  106  shows  the  ratio  of  iis  FiSD  values 
for  normal/tumor  tissues  at  each  measurement  wave¬ 
length.  A  RSD  ratio  of  1  implies  that  tumor  and  normal 
tissue  have  equivalent  structural  variation.  Interestingly,  in 
the  case  of  invasive  cancer,  the  normal /tumor  RSD  is 
consistently  greater  than  1  for  all  wavelengths.  This 
suggests  greater  scattering  spatial  heterogeneity  in  the 
well -differentiated  normal  tissue  than  in  the  lesion.  These 
observations  are  consistent  with  the  idea  that  disease 
progression  leads  to  loss  of  differentiation,  and  that 
differentiation  provides  the  fundamental  origin  for  spatially 
varying  scattering  signals. 

RSD  ratios  of  about  0.5  and  less  are  consistently 
observed  with  malignant  DCIS  and  benign  fibroadenomas. 
In  the  case  of  the  fibroadenoma,  this  is  likely  due  to 
enhanced  structural  variation  (versus  normal  tissue)  intro¬ 
duced  by  well-defined  zones  of  dense  collagen  fiber 
synthesis  and  glandular  shrinkage.  DCIS  appears  to  be 
midway  between  these  two  extremes,  possibly  due  to  the 
absence  of  invasion  and  general  preservation/stimulation  of 
extracellular  matrix.  These  data  further  confirm  the  sensitiv¬ 
ity  of  FDPM  to  the  cell  and  matrix  alterations  seen  in  Figure  7. 
However,  by  measuring  spatially  varying  properties,  the 
perturbative  effect  tumors  have  on  normal  tissue  architecture 
is  detected  with  enhanced  sensitivity. 


Conclusions 

FDPM  is  a  non-invasive  optical  technique  that  utilizes  NIR 
light  to  monitor  physiology  in  bulk  tissues.  Optical  proper-' 
ties  derived  from  FDPM  measurements  can  be  used  to 
construct  low- resolution  functional  maps  and,  conse¬ 
quently,  provide  a  relatively  low-cost  adjunct  to  many 
conventional  diagnostic  tools.  Substantial  work  remains  to 
accurately  describe  heterogeneous  tissues  that  vary  in 
geometry,  structure,  and  composition.  Despite  these 
limitations,  several  unique  applications  of  photon  migration 
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measurements  are  currently  underway,  particularly  in  solid 
tumor  characterization.  Preliminary  studies  show  that 
unique  breast  tumor  optical  signatures  are,  in  fact, 
detectable  in  human  subjects.  Changes  in  blood  flow, 
oxygen  consumption,  and  tissue  structure  (cellular  and 
matrix)  can  be  measured  quantitatively.  Most  importantly, 
these  signals  appear  to  reflect  fundamental  physiological 
processes  associated  with  malignancy,  disease  progres¬ 
sion,  age  and  menopausal  status.  Because  many  attributes 
related  to  the  origin,  progression,  and  therapeutic  response 
of  breast  cancer  are  common  to  other  types  of  tumors,  we 
expect  continued  optical  technology  development  will  have 
a  broad  impact  on  optical  diagnosis  and  therapy. 
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Abstract 


A  technique  for  measuring  broadband  near-infrared  absoprtion  spectra  of  tnrbid 
media  is  presented  using  a  combination  of  frequency-domain  (FD)  and  steady-state' 
(SS)  reflectance  methods.  Most  of  the  wavelength  coverage  is  provided  by  a  white- 
light  SS  measurement,  while  the  FD  data  are  acquired  at  a  few  selected  wavelengths. 
Coefficients  of  absorption  (fia)  and  reduced  scattering  (/r')  derived  from  the  FD 
data  are  Tised  to  intensity-calibrate  the  SS  measurements  and  to  estimate  at  all 
wavelengths  in  the  spectral  window  of  interest.  After  these  steps  are  performed, 
Ha  can  be  determined  by  comparing  the  SS  reflectance  values  to  the  predictions  of 
diffusion  theory,  wavelength  by  wavelength.  Absorption  spectra  of  a  turbid  phantom 
and  of  human  breast  tissue  in  vivo,  derived  using  the  combined  SSFD  technique, 
agree  well  with  expected  reference  values.  All  measurements  can  be  performed  at 
a  single  source-detector  separation  distance,  reducing  variations  in  sampling  volume 
that  exist  in  multi-distance  methods.  The  technique  uses  relatively  inexpensive  light 
sources  and  detectors  and  is  easily  implemented  on  an  existing  multi-wavelength  FD 
system. 

OCIS  codes;  170.1470,  Blood/tissue  constituent  monitoring;  170.4090,  modula¬ 
tion  techniqties;  170.5280,  photon  migration;  170.7050,  turbid  media. 
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1  Introduction 


Reflectance  spectroscopy  has  become  a  widely  used  technique  for  characterizing  turbid  me¬ 
dia,  including  many  medical  applications  to  human  tissue.  In  many  cases,  the  quantification 
of  chromophore  concentrations  is  desired,  and  this  requires  the  ability  to  separate  the  ef¬ 
fects  of  absorption  from  those  of  scattering.  Fundamentally,  the  coefficients  of  absorption, 
fia,  and  reduced  scattering,  /x'^,  can  be  determined  by  a  series  of  reflectance  measurements 
performed  in  one  of  three  domains,  namely  time^"^  (using  a  fast  pulse  of  light),  frequency^”^ 
(using  a  sinusoidally  modulated  source  of  light),  or  steady  state*"^^  (using  a  source  of  con¬ 
stant  intensity  but  multiple  detectors  at  different  distances).  Unsurprisingly,  these  three 
techniques  have  different  merits  and  limitations.  Spatially-resolved  steady-state  techniques 
are  relatively  inexpensive,  and  more  readily  suited  for  the  determination  of  and  /x(,  over 
large,  continuous  ranges  of  wavelengths  than  the  other  methods.  However,  the  steady-state 
approach  works  best  when  measurements  are  performed  using  a  combination  of  both  short 
(~  1  transport  mean  free  path)  and  long  (many  transport  mean  free  paths)  source-detector 
separations.^^  Ideally,  the  optical  properties  of  the  sample  should  not  vary  over  the  range  of 
different  volumes  probed  by  the  different  measurements.  The  larger  the  spread  of  distances 
probed,  the  more  likely  that  heterogeneities,  such  as  those  found  in  biological  tissue,  will 
distort  the  data  from  the  predictions  of  the  model.  One  approach  to  limiting  this  effect, 
given  that  the  shortest  separations  provide  great  stability  to  the  calculation  of  /x'^,  is  to 
use  relatively  short  source-detector  separations  (<10  mm).^^’^^  Since  the  mean  probing 
depth  scales  with  the  source-detecor  separation,  this  makes  such  measurements  sensitive 
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to  superficial  components  (to  depths  of  less  than  5  mm  for  typical  biological  tissues). 

Time-  and  frequencj'-domain  techniques  are  well  suited  for  deeper  investigations  (>1 
cm  for  biological  tissue).  Moreover,  they  can  be  performed  using  only  one  or  a  few  source- 
detector  separations,  which  makes  them  more  robust  in  studying  heterogeneous  samples. 
Because  such  techniques  require  sources  that  can  be  pulsed  or  modulated  rapidly,  covering 
a  large  wavelength  range  requires  a  tunable  laser  or  an  extensive  collection  of  laser  diodes, 
both  of  which  can  be  expensive,  difficult  to  maintain,  and  slow  to  cover  the  entire  spectrum. 
This  is  an  important  drawback,  since,  as  discussed  by  Hull  et  the  quantification  of 
chromophore  concentrations  can  be  significantly  affected  by  using  a  limited  number  of 
wavelengths. 

In  this  article,  we  suggest  a  way  to  use  steady-state  (SS)  and  frequency-domain  (FD)  re¬ 
flectance  measurements  in  tandem  to  obtain  broad  wavelength  coverage  with  increased  pen¬ 
etration  depth.  This  method  is  especially  promising  for  tissue  near-infrared  spectroscopy, 
such  as  breast  physiology  characterization.  For  such  applications,  the  method  proposed 
here  enables  rapid  data  acquisition,  deep  tissue  probing,  and  robust  resolution  of  the  con¬ 
tributions  from  the  four  major  NIR  tissue  absorbers:  oxy-  and  deoxy-hemoglobin,  water, 
and  fat.  The  central  innovations  are  a)  using  FD-derived  Ha  and  values  to  convert  the 
SS  measurements  into  units  of  absolute  reflectance,  and  b)  using  the  power-law  wavelength 
dependence  of  /r'  to  obtain  interpolated  and  extrapolated  values  at  non-laser  wavelengths. 
FD  measurements  are  performed  at  a  handful  of  diode  laser  wavelengths  spanning  the 
range  of  interest  (650-1000  nm),  while  the  SS  measurements  are  performed  continuously 
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across  the  entire  range.  Unlike  spatially-resolved  SS,  however,  only  a  single,  large  source- 
detector  separation  is  used,  preferably  the  same  one  as  for  the  FD  measurements.  The 
instrumentation  is  straightforward  and  particularly  easy  to  add  to  an  existing  FD  system. 

After  describing  our  method  of  combining  the  SS  and  FD  methods  (“SSFD”),  we  test 
the  SSFD  technique  by  using  it  to  measure  the  absorption  spectra  of  turbid  samples.  First, 
we  analyze  a  liquid  tissue  phantom  whose  absorption  spectrum  is  known  (via  direct  spec¬ 
trophotometry  of  its  absorbing  component  prior  to  mixing)  and  to  which  the  SSFD  result 
can  be  compared.  We  also  measure  locations  on  the  breast  of  two  human  female  volunteers, 
demonstrating  that  data  can  be  gathered  in  vivo  and  analyzed  to  provide  pertinent  phys¬ 
iological  parameters.  Differences  among  the  measured  breast  spectra  can  be  interpreted 
in  terms  of  different  relative  levels  of  water,  fat,  and  hemoglobin  present  in  the  explored 
tissue  volumes.  In  addition,  estimates  of  absolute  concentrations  are  comparable  to  those 
reported  in  other  recent  broadband  in  vivo  studies.  Lastly,  we  discuss  possible  reasons 
for  imperfect  spectral  fits,  and  we  compare  concentration  predictions  derived  from  SSFD 
analysis  to  those  derived  from  FD  data  alone. 

2  Experimental  methods 

2.1  Optical  measurements 

Figure  1  shows  the  experimental  arrangement  for  SSFD  measurements.  In  all  cases,  light 
is  delivered  via  optical  fiber  to  the  surface  of  the  sample  and  collected  at  some  distance  p 
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away.  For  the  liquid  phantom  measurements,  p  was  15.5  mm;  for  the  breast,  21.5  (FD)  and 
24  (SS)  mm  (the  slight  difference  was  due  to  instrumental  limitations;  a  future  instrument 
will  utilize  identical  distances).  In  FD  mode  (upper  box),  the  light  arrives  sequentially 
from  one  of  seven  amplitude-modulated  diode  lasers  (672,  800,  806,  852,  896,  913,  and  978 
nm,  all  with  output  powers  of  <  20  mW  at  the  sample)  and  is  detected  by  an  avalanche 
photodiode  unit  (Hamamatsu  C556P-56045-03)  that  amplifies  the  AC  component  of  the 
signal.  A  network  analyzer  (Hewlett  Packard  8753C)  delivers  251  modulation  frequencies 
between  100  and  700  MHz  and  measures  phase  and  modulation  amplitude  of  the  photon 
intensity  signal,  as  described  elsewhere.^*'®  In  SS  mode  (lower  box),  light  comes  from  a 
150W  halogen  lamp  (Fiber-Lite)  and  is  analyzed  via  a  fiber-coupled  spectrograph  (Ocean 
Optics  S2000)  with  a  linear  CCD  detector  between  525  and  1155  nm,  with  the  useful 
range  for  our  experiments  being  650-1000  nm.  The  spectrograph  records  a  total  df  2048 
points  (0.35  nm/pixel),  and  the  spectral  resolution  is  5  nm  (full  width  at  half  maximum). 
Light  is  delivered  to  the  sample  using  a  bundle  of  four  fibers  (bundle  diameter  600  pm) 
and  collected  using  a  single  fiber  of  diameter  1  mm.  The  spectrum  of  the  light  source 
is  measured  separately  by  inserting  the  source  and  detector  fibers  into  different  ports  of 
an  integrating  sphere  (Labsphere,  IS-040-SF).  Relative  reflectance  is  calculated  to  be  the 
sample  spectrum  divided  by  the  source  spectrum  (note  that  both  measurements  use  the 
same  delivery  fiber,  collection  fiber,  and  detector  apparaturs).  Total  acquisition  time  per 
sample  for  SSFD  measurements  is  on  the  order  of  40  s  (30  for  FD  and  10  for  SS).  Calculation 
of  Pa{^)  was  performed  according  to  the  methods  of  the  next  section  using  in-house  Matlab 
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(The  Math  Works,  Inc.)  code,  making  use  of  the  optimization  toolbox. 

2.2  Samples 

The  phantom  contained  17  mg  of  green  dye  (naphthol)  and  50  ml  of  an  aqueous  scattering 
suspension  (Intralipid-20%,  Pharmacia,  Inc.)  dissolved  in  950  ml  of  water.  Prior  to  the 
addition  of  Intralipid,  the  absorbance  spectrum  of  the  dye-water  solution  was  measured  in 
a  1  cm  plastic  cuvette  via  spectrophotometer  (Beckman  Instruments,  DU630)  with  water 
as  a  reference.  Optical  measurements  were  performed  as  indicated  in  Figure  1,  with  source 
and  detector  fibers  placed  at  the  surface  of  the  samples. 

In  vivo  measurements  were  performed  on  the  left  breasts  of  two  supine  female  volun¬ 
teers,  ages  37  and  21.  Data  were  gathered  from  two  regions  on  one  volunteer  (a  region  on 
the  areolar  border  and  a  region  of  “inner  breast”,  i.e.  close  to  the  middle  of  the  chest)  and 
one  region  on  the  other  volunteer  (inner  breast  only),  resulting  in  three  total  samples.  In 
this  modality,  the  source  light  was  again  delivered  by  optical  fiber  but  the  FD  detector  was 
placed  directly  against  the  tissue,  without  a  collection  fiber.  Fiber  and  detector  were  bun¬ 
dled  into  a  single  handheld  device  that  was  placed  gently  against  the  breast.  SS  reflectance 
was  measured  subsequently,  at  the  exact  same  location,  in  the  two-fiber  mode  described 
above.  All  procedures  were  approved  by  the  Institutional  Review  Board  of  the  University 
of  California,  Irvine  (study  95-563). 
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3  Background  Theory 


3.1  Diffusion  model 

Tlie  measured  reflectance  signal  R  (see  Figure  1)  can  be  predicted  theoretically  using 
the  diffusion  approximation  to  the  radiative  transfer  equation,  as  many  groups  have  dis¬ 
cussed/’ In  this  theoretical  framework,  the  reflectance  is  a  function  of  the  optical  prop¬ 
erties  of  the  medium,  defined  by  the  absorption  coefficient  ^a-,  reduced  scattering  coefficient 
/i(,  and  index  of  refraction  n.  The  diffusion  approximation  is  valid  for  large  source-detector 
separation  r  [r  >  10(/ra  +  and  high  (reduced)  albedo  >  0.95].  The 

SS  and  FD  cases  can  be  described  using  a  single  formalism  in  which  the  solution  for  the 
reflectance  is  built  from  the  Green’s  function  for  the  diffusion  equation,  i.e.  the  fluence 
[W/cm^j  due  to  an  isotropic  point  source  in  an  infinite,  homogeneous  medium.  This  func¬ 
tion  takes  the  form  exp  {—kr)/{Dr),  where  D  =  [3(/Za  +  the  complex  wavenumber 

k  —  k^eal  T  '^'kimagi  and 


with  u)  being  the  modulation  frequency  in  radians/sec.  The  steady-state  solution  is  simply 

the  limit  when  a;=0.  In  this  case,  the  solution  for  J?  is  a  real  number  {kreai  = 

and  kimag  —  0).  When  u  >  0,  R  becomes  a  complex  number  ylexp  {—icf)),  where  A  is  the 
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modulation  amplitude  and  (j)  is  the  phase  shift  relative  to  the  source.  As  explained  in  next 


section,  these  quantities  can  be  obtained  from  the  FD  measurement. 

In  treating  reflectance  problems,  the  sample/air  interface  is  modeled  using  an  ex¬ 
trapolated  boundary  condition, in  which  the  fluence  is  set  to  zero  at  a  distance  zi,  = 
2D{1  -f  Rcf[)/{1  —  above  the  sample  (/2cff  depends  upon  the  refractive  mismatch  and 
equals  0.493  for  tissue  of  n  =  1.4  and  air  of  n  =  1.0).  The  method  of  images  is  employed, 
with  an  isotropic  point  source  at  a  depth  Zo  =  (yUa  +  contributing  a  signal  Sr,  and 
a  negative  image  point  at  a  height  Zo  +  Zb  above  the  extrapolated  boundary  contributing 
Si-  This  gives  the  fluence  <I>  at  any  point  in  the  sample.  The  detected  signal  along  the 
boundary  is  then  written  as  a  combination  of  terms  proportional  to  the  fluence  and  to  its 
flux  normal  to  the  surface: 

R  =  Ci^  —  •  {—z),  (1) 

where  values  for  the  constants  Ci  and  C2  are  determined  by  the  refractive  index  mismatch 
between  the  two  media^^  (for  tissue  of  n  =  1.40  and  air  of  n  =  1,  these  constants  assume 
the  values  of  0.118  and  0.306,  respectively),  —z  is  a  unit  vector  pointing  normally  upwards 
out  of  the  sample,  and 


$  = 


P 

A-nD 


exp  {—krs)  exp  {—kvi) 


•  z  =  — 
47r 


rs  ri 

1  \  exp  {-kr^) 


+  {^o  +  2^6)  I  k  + 


1  \  exp  {—kri) 


,  (2) 


where  P  is  the  incident  power  and  Vg  (rj)  is  the  distance  from  the  source  (image)  to  the 
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detector. 


3.2  Frequency  domain  fitting 

Like  the  theory  of  the  previous  section,  our  frequency-domain  fitting  process,  which  provides 
Ha  and  h's  values  at  a  few  wavelengths,  has  been  described  at  length  elsewhere,^’ and  is 
provided  here  as  background. 

Each  FD  measurement  contains  instrumental  artifacts;  i.e. ,  the  measured  reflectance 
is  actually 

R^CoAexp[-i{(f)  +  (j)o)],  (3) 

where  Co  and  (j)o  are  sample-independent  instrumental  constants.  A  and  (j)  are  respectively 
the  modulation  amplitude  and  phase  of  the  frequency-domain  reflectance,  as  defined  in  the 
previous  section. 

For  calibration,  we  gather  FD  data  from  a  prepared  sample  whose  Ha  and  h's  values  are 
known  a  priori  from  a  set  of  two-distance  frequency-domain  measurements.^  Frequency- 
dependent  values  of  (po  and  Co  are  calculated  from  the  discrepancies  between  measured 
(Eq.  3)  and  predicted  (Eq.  1)  phase  and  modulation  amplitude,  thus  calibrating  our  future 
measurements. 

With  the  instrumental  constants  thus  determined,  the  502  data  points  per  sample  (251 
each  for  both  phase  and  amplitude  vs.  w)  depend  upon  two  unknowns:  pa  and  p's-  We 
select  the  best  pa  and  p'^  values  to  fit  the  predictions  of  Eq.  1  to  the  data  using  the  iterative, 
nonlinear,  least-squares  method  of  Levenberg  and  Marquardt,^^  alternating  between  fitting 
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the  j)hase  and  the  amplitude  until  overall  convergence  is  achieved/  As  noted  by  various 
authors, the  reflectance  in  Eq.  1  is  very  nearly  proportional  to  either  the  flux  or  the 
fluence  term  alone  (i.e. ,  the  ratio  of  the  two  is  independent  of  iia  and  /r'^)  for  the  large 
source-detector  separations  that  we  employ  (>  10  {/ia  +  We  therefore  use  only 

the  fluence  term  $  from  Eq.  1,  absorbing  the  additional  proportionality  factor  into  the 
instrumental  contstant  Co- 

4  Determination  of  broadband  fia  spectrum 

The  previous  section  provided  background  information  regarding  diffusion  theory  models 
and  a  method  for  extracting  /Xq  and  /x'  values  from  FD  measurements  at  specific  wave¬ 
lengths.  We  now  describe  how  FD  and  SS  measurements  can  be  combined  to  obtain 
quantitative,  broadband  /da  spectra. 

Our  goal  is  to  compute  /x^  at  each  wavelength,  given  the  measured  SS  reflectance. 
However,  single-distance  SS  reflectance  cannot  itself  provide  fia-  R  depends  upon  /x'^  as  well 
as  /ia,  so  one  measurement  cannot  provide  a  unique  determination  of  either  parameter.  In 
addition,  the  instrumental  constant  Co  is  not  known.  The  single  SS  measurement  at  each 
wavelength  therefore  needs  to  be  supplemented  with: 

1.  A  means  of  establishing  the  absolute  reflectance  intensity. 

2.  An  additional  piece  of  information  about  /Xa  or  (or  a  combination  of  the  two). 
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4.1  Use  of  FD  data  to  provide  additional  information 

The  FD  system  can  provide  both  pieces  of  necessary  information  at  each  SS  wavelength.  At 
first  glance  this  is  surprising,  as  the  FD  system  operates  at  only  seven  wavelengths,  while 
the  SS  system  covers  450  nm  with  continuous  5  nm  resolution.  It  would  seem  then  that 
the  needed  values  of  and  €„  could  be  determined  only  at  these  seven  wavelengths,  and 
strictly  speaking  this  is  indeed  the  case.  Significantly,  however,  the  wavelength  dependence 
of  both  fi'g  and  Co  are  smooth  and  predictable  in  shape,  as  will  be  discussed  below.  It  is 
this  fact  that  enables  us  to  use  the  discrete  FD  information  to  supplement  our  broadband 
SS  measurements  at  all  desired  wavelengths. 

Amplitude  calibration  of  reflectance  spectrum  The  FD  measurements  enable  cal¬ 
culation  of  the  instrumental  factor  Co  at  all  wavelengths  because  we  expect  no  wavelength 
dependence  at  all:  the  sample  and  source  spectra  are  measured  using  the  same  delivery  and 
collection  system  and  only  a  few  seconds  apart,  so  there  should  be  negligible  wavelength- 
dependent  artifacts  in  the  ratio  of  the  two.  As  such,  the  task  reduces  to  that  of  calculating 
Co  at  a  single  wavelength.  This  is  readily  done  at  any  of  the  FD  wavelengths,  as  //q  and 
fi'g  are  both  known,  and  absolute  reflectance  is  a  function  of  only  these  variables  (cf.  Eqs.  1 ' 
and  2  in  the  a;  =  0  limit).  Using  all  of  the  FD  data  to  increase  robustness,  we  can  calculate 
the  value  of  Co  that  scales  the  measured  SS  reflectances  to  match  the  predicted  reflectances 
as  closely  as  possible  (in  the  least- squares  sense).  Once  this  is  done,  the  scale  factor  Co  is 
known  for  the  entire  spectrum.  Examples  of  this  scaling  will  be  shown  below. 
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Determination  of  The  FD  values  of  /i'  also  allow  us  to  obtain  information  about 

/i^  across  the  entire  wavelength  range.  As  numerous  groups  have  observed, the  particle 
size  distribution  of  scatterers  (0.1-10/im)  in  many  biological  media  and  phantoms  tends 
to  have  smooth  wavelength  dependence  over  the  range  650-1000  nm,  well-described  by  a  - 
power  function  of  the  form 


IJ-sW  —  AX  ^ 


(the  parameter  A  here  having  no  relation  to  the  modulation  amplitude).  As  a  result, 
measuring  several  values  allows  us  to  fit  them  to  a  simple  function  of  wavelength  and 
obtain  good  estimates  of  at  all  other  wavelengths  needed. 


Iterative  solving  for  /ia(A)  spectrum  With  //^(A)  calculated  and  the  reflectance  cor¬ 
rectly  scaled,  the  equation  for  steady-state  reflectance  (Eq.  1)  contains  only  one  unknown, 
Ha-  Because  it  is  difficult  to  obtain  an  analytical  solution  for  jia  from  this  equation,  nu¬ 
merical  solving  methods  were  employed.  Proceeding  one  wavelength  at  a  time,  we  used  the 
MATLAB  nonlinear  zero-finding  function  fzero  to  choose  Ha  such  that  Rmeas  —  Rthy  it^a)  = 
0,  where  Rmeas  is  the  calibrated  SS  reflectance  measurement  and  Rthy  (l^a)  is  the  theoretical 
reflectance  predicted  by  Eqs.  1  and  2  for  a  given  trial  value  of  Solving  for  for  the  entire 
Ha  spectrum  required  about  10  seconds  on  a  personal  computer. 
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5  Results 


Figure  2  displays  the  SS  reflectance  measured  from  the  phantom  and  the  predicted  absolute 
reflectance  calculated  from  diffusion  theory  based  upon  the  FD  measurements  of  and  /r'^. 
Because  of  strong  absorption  at  wavelengths  above  950  nm,  the  FD  measurement  at  978 
nm  was  very  noisy  at  short  integration  times,  which  in  turn  made  calculation  of  SS  data 
at  wavelengths  above  950  unstable.  Data  are  therefore  presented  between  650-950  nm,  for 
the  phantom  only.  As  anticipated,  all  of  the  measurements  differ  from  the  predictions  by 
essentially  the  same  scale  factor.  The  entire  SS  spectrum  is  thus  converted  to  absolute 
reflectance  units.  Similar  results  were  observed  for  the  three  breast  samples,  as  shown  in 
Figure  3.  In  the  order  presented,  sample  1  is  the  inner  breast  region  of  the  37-year-old 
volunteer,  and  samples  2  and  3  are  respectively  from  the  areolar  border  and  inner  breast  ■ 
of  the  21-year-old  volunteer. 

Figure  4  shows  the  least-squares  power-law  fit  to  the  values  measured  by  FD  on  the 
phantom.  As  expected,  the  wavelength  dependence  is  fairly  smooth  and  easily  described 
by  the  fit.  Corresponding  fits  for  the  three  breast  samples  appear  in  Figure  5.  Note  that 
the  values  only  change  by  a  factor  of  2  over  the  entire  range,  so  the  wavelength  dependence 
is  weak  in  addition  to  being  predictable. 

Figure  6  summarizes  the  absorption  values  calculated  for  the  dye-Intralipid  phantom. 
The  thick  line  shows  the  Ha  values  that  were  found  to  make  Eq.  1  reproduce  the  measured 
reflectance  data  as  closely  as  possible,  using  the  fzero  algorithm  as  previously  described. 
The  thin  line  represents  a  linear  least-squares  fit  to  the  thick  line  based  upon  refetence 
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spectra  of  pure  water  (from  Kou  et  and  pure  dye  (measured  by  spectrophotometer 

as  described  above).  As  the  figure  shows,  the  two-component  fit  accurately  models  the 
measured  spectrum  across  the  entire  spectral  range.  The  least-squares  fitting  coefficient 
for  the  dye  spectrum  is  1.00.  Also  shown  are  the  values  of  jig.  derived  by  FD  analysis  alone, 
demonstrating  an  agreement  between  the  FD  and  SSFD  methods  at  these  wavelengths,  as 
would  be  expected  because  the  FD  data  have  been  used  to  calibrate  the  SS  data. 

Corresponding  plots  showing  absorption  spectra  of  the  three  breast  samples  appear  in 
Figures  7  through  9.  As  before,  the  thick  line  is  the  experimental  data,  the  thin  line  is  a 
full-spectrum  fit,  and  the  circles  are  FD  values.  In  these  cases,  however,  the  fit  is  built 
from  published  spectra  (see  Figure  10)  of  oxy-  and  deoxy-hemoglobin,^^  water,^®  and  fat,^® 
which  are  commonly  regarded  as  the  four  major  absorbers  in  breast  tissue  in  the  650-1000 
nm  regime.  Coefficients  from  the  fit  thus  provide  estimates  of  these  components’  con¬ 
centrations.  The  estimates  are  useful  for  comparing  one  sample  to  another;  no  attempt  at 
independent  concentration  validation  has  been  performed  (however,  they  are  approximately 
correct  in  magnitude  and,  as  noted  in  the  next  section,  the  values  are  consistent  with  others 
from  the  recent  literature^’ ^®).  For  comparison  purposes,  the  additional  dotted  line  shows 
the  spectral  fit  obtained  when  only  the  FD  wavelengths  are  weighted,  as  in  a  customary 
FD-only  experiment.  The  two  spectral  reconstructions  are  clearly  different,  which  leads 
to  different  concentration  estimates.  All  of  the  concentration  predictions,  both  SSFD  and 
FD,  are  listed  in  Table  1,  along  with  a  final  column  indicating  the  percent  deviation  of 
the  FD  value  from  the  SSFD  value.  The  hemoglobin  results  are  reported  in  terms  of  total 
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hemoglobin  (oxy  plus  deoxy)  and  oxygen  saturation  (oxy  divided  by  total). 

6  Discussion 

In  order  to  extract  fia  from  a  single  SS  reflectance  spectrum,  one  has  to  convert  the  data  into 
absolute  units  {i.e. ,  the  instrumental  factor  Co  must  be  removed)  and  determine  values 
at  all  wavelengths.  We  have  demonstrated  that  FD  measurements  at  a  few  wavelengths  can 
accomplish  these  goals  across  the  entire  650-1000  nm  spectrum.  This  somewhat  surprising 
result  arises  because  follows  a  power  law,  while  Co  has  no  wavelength  dependence 

at  all.  Two  innovations  that  combine  FD  and  SS  methods  are  thus  available.  First,  by 
scaling  SS  reflectance  to  match  FD  predictions  of  absolute  reflectance,  the  entire  reflectance 
spectrum  is  automatically  calibrated.  Secondly,  by  fitting  the  FD  values  to  a  wavelength- 
dependent  function,  the  extra  information  needed  to  extract  i^a  from  diffusion  theory  (Eqs.  1 , 
and  2)  is  obtained.  We  note  that  the  FD-derived  Ha  and  values  could  just  as  readily  be 
supplied  by  time-domain  measurements  at  a  single  source-detector  separation. 

The  interpolation  of  the  reduced  scattering  spectra  by  using  a  fit  is  an  important  point 
of  this  technique.  The  lower  two  plots  in  Figure  5  show  excellent  agreement  with  the 
power  law.  The  top  plot  shows  reasonably  good  agreement,  but  small  discrepancies  can 
nevertheless  be  observed  (on  the  order  of  5-8%).  This  behavior  is  likely  not  to  be  an 
actual  divergence  from  the  power  law,  but  is  probably  due  to  a  coupling  effect  between  the 
scattering  and  absorption  quantification  in  the  FD  fit.  Such  a  coupling  appears  when  the 
diffusion  model’s  assumptions  are  not  sufficiently  fulfilled.  The  tissue  heterogeneity,  i.e.  the 
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layered  superficial  structure  and  the  mix  of  deeper  fat  lobules  and  fibrous  tissue,  could  be 
a  possible  source  of  deviation  from  the  model.  Interestingly,  the  highest  fat  concentration 
was  measured  for  this  breast  location  (cf.  Table  1).  As  described  recently  by  Doornbos  et 
using  a  power  law  fit  to  calculate  the  spectrum  can  in  fact  provide  an  advantage 
in  computing  Ha-  Indeed,  by  smoothing  the  spectrum,  it  reduces  the  spurious  coupling 
between  scattering  and  absorption  caused  by  inaccurate  modeling. 

When  compared  to  discrete-wavelength  FD  measurements,  the  access  to  a  continuous, 
absorption  spectrum  offers  two  important  advantages:  better  chromophore  identification 
and  improved  concentration  quantification.  The  ability  of  the  method  to  recover  the  con¬ 
tinuous  shape  of  the  true  absorption  spectrum  is  evidenced  by  the  phantom  measurement. ' 
Figure  6  shows  that  the  experimental  spectrum  is  accurately  fitted  using  known  spectra 
of  dye  and  water.  The  fitting  coefficient  of  1.00  for  the  dye  additionally  demonstrates  the 
accurate  recovery  of  chromophore  concentrations. 

The  breast  spectra  (Figures  7  through  9)  further  illustrate  the  usefulness  of  the  SSFD 
technique  for  clinical  investigations.  As  expected,  the  fit  using  oxy-and  deoxy-hemoglobin, 
fat,  and  water  accounts  for  most  of  the  absorption  in  the  650-1000  nm  in  breast  tissues.  The 
quality  of  the  fit  is  especially  striking  for  wavelengths  longer  than  800  nm.  In  particular, 
the  fat  spectrum  used  in  the  fit  seems  accurate,  matching  the  main  peak  at  928  nm  and 
the  shoulder  at  895  nm.  These  results  are  surprisingly  good  considering  that  the  pure 
fat  spectrum  was  measured  on  soybean  oil.^*  These  measurements  stress  the  important 
contribution  of  fat  to  tissue  absorption,  as  this  chromophore  has  been  sometimes  neglected 
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in  previous  work.®’^'’  The  water  peak  at  976  nm  is  also  well  reproduced  in  the  measurements. 
We  observe  that  the  highest  water  content  was  measured  for  the  areolar  border  sample 
(sample  7^2).  Interestingly,  we  note  that  for  both  dye  and  breast  tissue,  data  fits  to  the 
water  spectrum  of  Kou  et  al.  are  superior  to  Hale  and  Querry,^®  particularly  in  the  920-960 
nm  regime  where  their  water  absorption  spectrum  increases  sharply. 

In  the  650-800  nm  region,  the  major  features  of  tissue  absorption  are  clearly  due  to 
oxy-  and  deoxy-hemoglobin.  The  deoxy-hemoglobin  peak  at  760  nm  is  distinctly  visible  in 
the  experimental  data.  Nevertheless,  we  observe  small  but  consistent  spectral  differences 
on  the  order  of  0.001  mm~^  between  the  fit  and  measurement  in  this  regime  (see  Figures  7 
through  9).  These  discrepancies,  revealed  in  the  full  SSFD  spectrum,  are  not  evident  when 
fitting  the  FD  data  alone.  Two  reasons  for  the  imperfect  fits  can  be  suggested.  First,  the 
oxy-and  deoxy-hemoglobin  absorption  spectra  we  used  could  be  slightly  incorrect.  Indeed, 
small  variations  are  found  between  various  published  spectra,  and  changes  could  be  also 
expected  between  in  vitro  and  in  vivo  values.  Second,  other  “background”  chromophores 
should  probably  be  taken  into  account.  For  example,  the  tails  of  several  protein  absorption 
bands,  not  included  in  the  modeling,  might  contribute  in  this  wavelength  region.  A  more 
extensive  spectral  library  is  therefore  desirable  for  future  studies.  Alternatively,  one  could 
model  the  background  empirically,  either  with  predetermined  mathematical  functions  of 
wavelength  {e.g.  Hull  et  or  via  a  principal-component  analysis  of  several  background 
residuals  once  many  samples  have  been  studied.^® 

For  the  sake  of  comparison,  we  used  the  same  least-square  fitting  algorithm  to  calculate 
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tlie  spectral  fits  and  concentrations  twice,  once  using  the  full  spectrum  and  once  using 
only  the  FD  wavelengths,  even  though  this  second  procedure  does  not  fully  exploit  the 
potential  of  the  SSFD  method.  Interesting  differences  are  found  between  the  FD  and 
SSFD  calculations.  Figure  7  and  Table  1  show  that  FD  significantly  underestimates  the 
fat  concentration  in  the  first  breast  sample:  0.75  g/cm^  using  SSFD  versus  0.44  g/cm^ 
using  FD.  The  reason  for  this  discrepancy  lies  in  the  fact  that  no  laser  wavelength  is  close 
enough  to  the  fat  peak  (at  present,  commercial  laser  diodes  at  this  particular  wavelength 
are  uncommon).  Thus,  with  no  weighting  in  this  region,  the  FD  fit  tolerates  large  fitting 
errors  near  the  absorption  maximum  of  fat  (as  seen  in  Figure  7),  leading  to  large  errors  in 
the  fat  concentration.  Such  a  discrepancy  is  naturally  enhanced  for  this  sample,  where  the 
fat  concentration  is  highest.  This  example  illustrates  the  shortcomings  of  using  a  limited 
number  of  sources  (be.  wavelengths),  especially  when  a  significant  chromophore  lies  in  an 
undersampled  specral  region.  A  similar  error  is  seen  in  the  fitting  of  the  areolar  breast 
location  (Figure  8),  where  this  time  it  is  the  water  peak  that  is  poorly  addressed  by  the 
FD  analysis. 

The  total  hemoglobin  determination  is  more  consistent  between  the  FD  and  SSFD 
calculations,  due  to  the  smoother  spectral  features  of  these  chromophores  and  better  diode 
coverage  of  the  650-850  nm  regime.  However,  the  oxygenation  determination  by  FD  is 
systematically  higher.  As  discussed  before,  consistent  small  differences  occur  between  the 
fit  and  the  experiments  in  the  650-800  nm  region  due  to  some  incompleteness  in  our  library 
of  fitting  lineshapes.  Fitting  the  curve  with  a  limited  number  of  wavelengths  naturally 
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exacerbates  the  incompleteness  and  produces  a  bias  in  the  results. 

As  mentioned  above,  Cubeddu  et  al?  obtained  similar  optical  properties  and  physiolog¬ 
ical  parameters  for  premenopausal  breast,  working  in  the  time  domain.  Conceptually,  the 
two  experiments  exploit  the  same  optical  phenomena  and  are  essentially  complementary 
techniques  in  different  domains.  A  technical  advantage  of  the  SSFD  measurement  is  the 
use  of  a  white  light  source  rather  than  a  tunable  laser.  As  such,  the  SS  measurement  cov¬ 
ers  all  wavelengths  simultaneously,  whereas  the  laser(s)  must  be  tuned  separately  to  each 
wavelength.^  Additionally,  the  heart  of  the  SSFD  system  (a  few  laser  diodes,  a  frequency 
generator,  an  avalanche  photodiode,  a  network  analyzer  or  lock-in  amplifier,  and  a  SS  re¬ 
flectance  system)  is  inexpensive  and  easy  to  maintain  compared  to  a  tunable  laser  system 
and  a  single-photon-counting  detection  apparatus. 

As  already  stated,  the  method  presented  here  has  the  advantage  of  being  compatible 
with  a  single,  large  source-detector  separation  (:^10  transport  mean  free  paths).  In  contrast 
to  spatially-resolved  methods,  it  is  well  suited  for  interrogating  deep  strucures  in  relatively 
heterogeneous  samples.  As  numerous  groups  have  shown^^’®^’^^  the  layered  structure  of 
tissue  affects  reflectance  differently  at  different  source-detector  separations,  raising  doubts 
about  the  applicability  of  spatially-resolved  techniques  that  assume  sample  homogeneity 
over  a  large  range.  While  variations  in  the  FD  modulation  frequency  w  do  change  the 
optically-explored  tissue  volume,  these  effects  are  modest  in  comparison  to  changes  in  p, 
particularly  at  high  absorption. Consequently,  the  assumption  of  homogeneity  is  less 
extreme  for  the  single-distance  measurements  proposed  here.  The  essence  is  that  all  mea- 
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surements  at  one  wavelength  explore  a  more  consistent  volume;  this  does  not  mean  that 
different  wavelengths  will  explore  the  same  volume.  Indeed,  less-absorbed  wavelengths  will 
explore  larger  regions  than  wavelengths  that  attentuate  more  rapidly. 

7  Conclusion 

A  combination  of  SS  and  FD  reflectance  measurements  has  been  described  for  absorption 
spectroscopy  of  turbid  media,  featuring  beneficial  aspects  of  both  techniques.  As  with  SS, 
the  wavelength  coverage  is  continuous,  detecting  absorption  features  that  may  elude  the  dis¬ 
crete  wavelengths  chosen  for  FD.  The  prediction  of  constituent  concentrations,  for  instance 
in  breast  tissue,  is  substantially  improved  using  full-spectrum  absorption  data  rather  than 
a  handful  of  wavelengths.  As  with  FD,  however,  only  a  single  source-detector  separation 
is  required,  making  the  technique  more  amenable  to  reporting  volume- averaged  values  for 
heterogeneous  samples.  In  addition,  the  source-detector  separation  can  be  large,  allowing 
for  centimeter-scale  mean  probing  depths  that  cannot  be  achieved  with  spatially-resolved 
SS  techniques.  This  advantage  for  deeply-probing  studies  is  significant  for  many  clinical 
purposes.  An  application  to  breast  analysis  has  been  demonstrated,  with  quantitative  in 
vivo  spectra  of  human  breast  obtained  rapidly  (<1  min.).  The  total  hemoglobin  content, 
oxygen  saturation,  water,  and  fat  content  of  the  breast  samples  have  been  calculated  from 
the  spectra,  and  failures  of  FD-only  fitting  have  been  highlighted.  The  technique  is  rela¬ 
tively  inexpensive  and  could  prove  valuable  for  improving  accuracy  in  the  development  of 
quantitative  photon  migration  for  clinical  instruments. 
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predictions  for  the  third  breast  measurement  (inner  breast,  21-year-old 
volunteer).  Lines  and  circles  have  the  same  assignments  as  described  in 
Figure  7.  Physiological  parameters:  total  hemoglobin  concentration,  19  /zM; 
oxygen  saturation,  72%;  water,  28  g/cm^;  and  fat,  0.56  g/cm^.  See  Table  1 

for  more  details .  39 

10  Absorption  spectra  of  major  absorbers  in  breast  tissue.  Upper  left:  oxyhe¬ 
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Figure  1:  Configuration  of  light  sources,  optical  fibers,  and  detectors. 
Key:  APD,  avalanche  photodiode.  The  two  dotted  rectangles  denote 
components  belonging  to  the  FD  and  SS  systems,  respectively. 
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Figure  2:  SS  reflectance  spectrum  acquired  from  dye-Intralipid  phantom, 
scaled  to  fit  the  discrete  reflectance  values  (circles)  predicted  by  Eqs.  1 
and  2  using  FD-derived  values  of  and  /ia.  This  scaling  causes  the 
entire  spectrum  to  be  converted  into  absolute  reflectance  units;  see  text 
for  discussion.  Error  in  FD  reflectance  is  estimated  at  ±3  x  10“^/mm; 
in  SS  reflectance,  ±1  x  10“^/mm. 
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Figure  3:  SS  reflectance  spectra  acquired  from  three  locations  on  breast 
tissue  of  female  volunteers.  The  spectra  have  been  scaled  to  fit  reflectance 
values  (circles)  calculated  from  FD  data  in  the  same  manner  as  Figure  2. 
Errors  are  same  as  previous  figure. 
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Figure  4:  Determination  of  broadband  /r'  spectrum  for  the  dye-Intralipid 
phantom.  Circles  are  the  discrete  values  measured  by  the  FD  tech¬ 
nique,  while  the  solid  line  is  the  best  power-law  fit  to  Eq.  4.  Fitting 
parameters:  A  =  2200,  B  =  -0.82.  Error  bars  on  are  3%.^ 
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Figure  5:  Determination  of  broadband  yit'  spectra  for  the  three  breast 
samples.  Circles  are  the  discrete  /x'  values  measured  by  the  FD  tech¬ 
nique,  while  the  solid  lines  are  the  best  power-law  fit  to  Eq.  4.  Fitting 
parameters  (A,B)  are  (240,-0.86),  (2700,-12),  and  (250,-0.85),  respec¬ 
tively.  Error  bars  on  are  3%.^ 
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Figure  6:  Comparison  of  Ha  values  generated  by  FD  (dots)  and  by  SSFD 
(thick  line)  for  the  dye-Intralipid  phantom.  Also  shown  is  the  best  fit 
(thin  line)  to  the  SSFD  data  using  spectra  of  naphthol  (measured  on 
a  spectrophotometer)  and  water  (from  Kou  et  Circles  are  the 

discrete  pia  values  from  FD  alone.  Error  bars  are  ±0.0005  mm“^  or  5%, 
whichever  is  larger. 
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Figure  7;  iia  predictions  for  the  first  breast  measurement  (inner  breast, 
37- year-old  volunteer).  The  thick  line  is  the  SSFD  data;  the  thin 
line  is  the  least-squares  fit  using  a  superposition  of  Hb,  Hb02,  wa¬ 
ter,  and  fat  spectra;  the  circles  are  the  FD  values  (error  bars  are 
±0.0005  mm“^  or  5%,  whichever  is  larger);  the  dotted  line  is  the  least- 
squares  fit  when  only  the  FD  values  are  weighted.  Physiological  param¬ 
eters  from  SSFD  spectral  fit:  total  hemoglobin  concentration,  22  /rM; 
oxygen  saturation,  73%;  water,  15  g/cm^;  and  fat,  0.75  g/cm®.  See  Ta¬ 
ble  1  for  more  details.  Note  that  the  fat  peak  between  900  and  950  nm 
is  significantly  underfit  by  the  FD  calculation. 
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Figure  8:  fia  predictions  for  the  second  breast  measurement  (areolar  bor¬ 
der,  21-year-old  volunteer).  Lines  and  circles  have  the  same  assignments 
as  described  in  Figure  7.  Physiological  parameters:  total  hemoglobin 
concentration,  30  /xM;  oxygen  saturation,  70%;  water,  0.51  g/cm^;  and 
fat,  0.42  g/cm^.  See  Table  1  for  more  details. 
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wavelength  (nm) 

Figure  9:  predictions  for  the  third  breast  measurement  (inner  breast, 

21-year-old  volunteer).  Lines  and  circles  have  the  same  assignments  as 
described  in  Figure  7.  Physiological  parameters:  total  hemoglobin  con¬ 
centration,  19  pM;  oxygen  saturation,  72%;  water,  28  g/cm^;  and  fat, 
0.56  g/cm^.  See  Table  1  for  more  details. 
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Figure  10;  Absorption  spectra  of  major  absorbers  in  breast  tissue.  Up¬ 
per  left:  oxyhemoglobin,  1  /rM;  upper  right,  deoxyhemoglobin,  1  /zM; 
lower  left,  water,  1  g/cm^;  lower  right,  fat  (soybean  oil),  0.9  g/cm®.  The 
hemoglobin  spectra  are  from  Wray  et  al.^^  the  water  is  from  Kou  et 
and  the  fat  is  from  the  doctoral  thesis  of  Eker.^* 


40 


List  of  Tables 


1  Results  of  physiological  predictions  for  breast  samples  1-3.  The  final  col¬ 
umn  gives  the  %  difference  between  the  FD  and  SSFD  values,  defined  as 
100x(FD-  SSFD) /SSFD . 


component 

SSFD 

FD  only 

%  difference 

total  Hb  (/xM) 

22 

25 

13 

breast  1 

O2  saturation  (%) 

73 

78 

7 

water  (g/cm^) 

0.15 

0.14 

-7 

fat  (g/cm^) 

0.75 

0.44 

-41 

total  Hb  (/xM) 

30 

30 

0 

breast  2 

O2  saturation  (%) 

70 

79 

13 

water  (g/cm^) 

0.51 

0.40 

-22 

fat  (g/cm^) 

0.42 

0.39 

-7 

total  Hb  (^M) 

19 

19 

0 

breast  3 

O2  saturation  (%) 

72 

81 

13 

water  (g/cm^) 

0.28 

0.27 

-4 

fat  (g/cm^) 

0.56 

0.54 

-4 

Table  1:  Results  of  physiological  predictions  for  breast  samples  1-3.  The 
final  column  gives  the  %  difference  between  the  FD  and  SSFD  values, 
defined  as  100x(FD-  SSFD)/SSFD. 
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ABSTRACT 


Purpose:  Near  infrared  (NIR)  diffuse  optical  spectroscopy  and  imaging  may  enhance 

existing  technologies  for  breast  cancer  screening,  diagnosis,  and  treatment.  NIR  techniques  are 
based  on  sensitive,  quantitative  measurements  of  functional  contrast  between  healthy  apd 
diseased  tissue.  In  this  study  we  quantified  the  origins  of  this  contrast  in  healthy  breast. 

Materials  and  Methods:  Non-invasive  NIR  measurements  were  performed  on  the  breasts  of 

28  healthy  women.  Subjects  included  pre-  and  post-menopausal  women  between  the  ages  of  18 
and  64.  A  diffusive  model  of  light  transport  quantified  oxygenated  and  deoxygenated 
hemoglobin,  water,  and  lipid  by  their  absorption  signatures.  Cellular  density,  fat,  and  collagen 
content  were  inferred  from  measured  light-scattering  spectra. 

Results:  Substantial  quantitative  differences  were  observed  in  the  optical  properties  of 

breast  as  a  function  of  age.  Changes  were  consistent  with  long-term  hormone-dependent 
transformations  that  occur  in  breast.  Instrument  response  was  not  adversely  affected  by  subject 
age  or  menopausal  status. 

Conclusions:  These  measurements  provide  new  insight  into  endogenous  optical  absorption  and 
scattering  contrast  mechanisms  and  have  implications  for  the  development  of  optical 
mammography.  NIR  spectroscopy  yields  quantitative  functional  information  that  cannot  be 
obtained  with  other  non-invasive  radiological  techniques.  The  prototype  bedside-capable 
instrument  is  comparable  in  cost  to  commercialized  ultrasound. 


Key  Words:  breast,  breast  radiography,  blood,  fat 
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INTRODUCTION 


A.  Scope  of  This  Work 

This  paper  describes  the  application  of  a  new  class  of  near-infrared  (NIR)  spectroscopy  tools  for 
non-invasive,  quantitative  characterization  of  breast  tissue  physiology.  Diagnostic  methods 
currently  in  use,  such  as  mammography,  magnetic  resonance  imaging,  and  ultrasound  offer 
excellent  anatomical  lesion-detection  capabilities,  but  are  generally  unable  to  provide 
quantitative  information  regarding  tissue  function  and  composition  [1].  Positron  emission 
tomography  shows  great  promise  in  evaluating  the  metabolic  demands  of  tissue,  but  requires 
exogenous  radionuclides  and  is  insensitive  to  tissue  hemodynamics.  NIR  optical  spectroscopy  is 
intrinsically  sensitive  to  the  principal  components  of  breast:  blood,  water,  and  adipose,  as  well  as 
epithelial  and  connective  tissues.  Preliminary  studies  suggest  that  the  fractional  contribution  of 
each  to  NIR  signals  depends  strongly  on  factors  such  as  age,  menopausal  status,  and  the 
progression  of  disease.  Thus,  NIR  optical  spectroscopy  provides  an  opportunity  for  revealing 
physiologieal  information  that  is  unobtainable  by  other  non-invasive  techniques. 

As  optical  diagnostic  methods  move  rapidly  towards  the  clinic,  there  is  a  growing  need  to 
identify  and  quantify  sources  of  endogenous  contrast  in  breast  [2-11].  This  paper  presents  the 
results  of  quantitative  NIR  spectroscopic  studies  in  a  population  of  28  subjects  ranging  in  age 
from  1 8  to  64.  We  provide  evidence  that  NIR  tissue  spectroscopy  is  inherently  sensitive  to  long¬ 
term  age  and  hormone-dependent  dynamics.  Similar  studies  have  been  initiated  [12-14]  but  on 
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only  a  few  subjects  or  with  limited  sensitivity  to  physiological  parameters.  These  data  are 
expected  to  play  an  important  role  in  characterizing  the  sensitivity  of  optical  spectroscopy  and 
imaging  for  applications  such  as  tumor  detection,  therapeutic  monitoring,  and  risk  assessment. 


B.  Near-Infrared  Tissue  Spectroscopy 

NIR  photons  (600-1000  nm)  are  non-ionizing  and  weakly  absorbed  by  tissue.  Many  have 
exploited  their  enhanced  penetrance  relative  to  that  of  visible  or  ultra-violet  radiation  to  explore 
non-invasively  into  brain,  breast,  and  muscle  [1].  In  breast  tissue,  the  principal  NIR’ absorbers 
are  assumed  to  be  reduced  hemoglobin  (Hb-R),  oxygenated  hemoglobin  (Hb-Oa),  water  (H2O), 
and  lipids  [15].  NIR  tissue  spectroscopy  offers  a  safe  way  to  quantify  these  components  and  to 
view  unique  functional  information  with  low-cost  “point-of-care”  devices.  Multiple  light 
scattering  complicates  quantitative  measurements  of  light  absorption  in  the  NIR.  The  process  of 
using  frequency-domain  photon  migration  (FDPM)  techniques  to  separate  the  effects  of 
absorption  from  scattering  has  been  described  in  detail  [16,  17].  FDPM  enables  the  use  of 
quantitative  spectroscopic  analysis  tools  to  determine  the  composition  and  structure  of  tissue 
non-invasively  [18]. 

The  application  of  FDPM  techniques  to  NIR  tissue  spectroscopy  differs  substantially  frorh 
conventional  bulk  NIR  methods,  such  as  diaphanography  [19].  FDPM  employs  computational 
models  that  allow  the  quantitative  separation  of  the  effects  of  absorption  and  scattering.  Thus, 
within  the  limits  of  applicability  of  the  light-transport  models  to  the  breast,  all  parameters 
measured  in  this  study  are  absolute  values. 
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MATERIALS  AND  METHODS 


A.  Patient  Selection 

All  volunteers  enrolled  in  this  study  competently  provided  informed  consent  for  participating  in 
one  of  two  trials  (#95-563  and  #99-2183)  under  the  guidelines  of  an  Institutional  Review  Board. 
The  28  volunteer  ages  ranged  between  1 8  and  64.  Fifteen  of  the  volunteers  were  pre-menopausal 
(average  age  28  +!-  9).  Seven  were  post-menopausal  (average  age  56  +/-  2).  The  remaining  six 
(average  age  56  +/-  5)  were  women  taking  some  form  of  hormone  replacement  therapy  (HRT); 
three  of  these  six  women  classified  themselves  as  peri -menopausal.  None  of  the  women  had  any 
known  cancerous  lesions  or  other  known  forms  of  breast  disease. 


B.  Measurement  Technique 

Each  volunteer  rested  in  a  supine  position  during  the  measurement.  The  instrument  probe,  which 
is  slightly  larger  than  an  ultrasound  probe  (Figure  1),  was  the  only  item  placed  in  contact  with 
the  volunteer.  This  probe  consisted  of  a  sealed  avalanche  photodiode  (APD)  module  and  an 
optical  fiber  placed  22  mm  away  from  the  APD.  All  measurements  were  performed  in  a 
reflection-style  geometry.  The  probe  was  placed  on  the  breast  with  minimal  pressure  using  only 
the  force  of  gravity;  no  compression  was  used.  In  this  configuration,  we  estimate  that  the  light 
sampled  approximately  1  cm  below  the  skin.  We  report  data  only  on  measurements  performed 
in  the  center  of  the  left  upper  outer  quadrant  (measurements  on  the  right  upper  quadrant  yielded 
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similar  results).  Error  bars  in  Figures  and  Tables  represent  the  standard  deviation  of  repeated 
measurements. 


C.  Instrumentation 

A  sehematic  of  our  1  GHz,  portable  FDPM  device  is  shown  in  Figure  1 .  The  specific  details  of 
this  FDPM  instrument  have  been  descrihed  in  detail  elsewhere  [20],  hut  the  relevant  technical 
information  is  mentioned  here.  The  instrument  employs  multiple  diode  lasers  that  provide 
visible  and  NIR  light  at  seven  wavelengths  (672,  800,  806,  852,  896,  913,  and  978  nm).  A  hand¬ 
held  probe  has  been  designed  to  house  an  APD  that  records  the  modulated  diffuse  light  signals 
after  propagation  through  the  tissue.  This  probe  has  a  plastic  attachment  on  the  casing  to 
position  a  source  optical  fiber  a  fixed  distance  fi'om  the  APD.  A  100-pm-diameter  graded  index 
optical  fiber  positioned  22  mm  from  the  APD  detector  was  used  to  deliver  the  diode  laser  outputs 
to  the  tissue  surface.  A  network  analyzer  measured  the  phase  and  amplitude  of  the  detected 
electronic  signal  from  the  APD.  A  steady-state  current  source  was  mixed  with  RF  power 
provided  by  the  network  anal5^er  in  a  bias  network.  This  bias  network  serially  distributed  power 
to  each  laser  diode  and  produced  intensity-modulated  light.  An  optical  switch  delivered  light 
serially  from  each  diode  to  the  tissue  via  the  single  optical  fiber  described  in  the  handheld  probe. 
The  optical  power  launched  into  the  tissue  ranged  from  5-25  mW  for  each  wavelength. 

A  sweep  over  all  seven  wavelengths  ranged  from  approximately  35  to  60  s.  The  system  acquired 
data  in  less  than  3  s  per  wavelength,  with  a  2  s  delay  between  each  laser  diode  in  the  system 
because  of  switching  considerations.  The  system  was  wheeled  into  a  medical  clinic  for  the  each 


8 


measurement.  Instrumental  artifacts  were  removed  by  calibrating  on  a  tissue-simulating 
phantom  with  known  absorption  and  scattering  properties. 


D.  Measured  Parameters  and  Data  Analysis 

The  amplitude  and  phase  of  a  NIR  diffusive  light  wave  demodulates  and  retards,  respectively,  as 
the  wave  propagates  through  multiple-scattering  media  such  as  tissue.  The  real  and  imaginary 
parts  of  this  diffuse  wave  were  fit  simultaneously  to  a  light-diffiision  model  (the  Pi 
Approximation  to  the  Boltzmann  Transport  Equation)  [21,22]  by  minimizing  the  chi-squared 
( X  )  surface  with  a  Marquardt-Levenberg  algorithm.  This  fit  determined  the  absolute  optical 
absorption  coefficient,  ]ia  and  the  absolute  optical  reduced  scattering  coefficient,  at  each 
wavelength.  Typical  jia,  and  uncertainties,  determined  from  the  distribution  of  the  fits,  are 
1%  to  3%  of  their  mean  values.  In  actual  tissue  measurements,  precision  error  is  less  than  10% 
for  Pa  and  5%  for  p^'. 

When  the  optical  properties  p^'  and  p^  are  recovered  for  the  seven  wavelengths,  the  spectral- 
dependence  of  the  absorption  may  be  combined  with  known  values  of  molecular  extinction 
coefficients  to  calculate  physiologically-relevant  parameters.  We  assume  that  bulk  breast  tissue 
is  composed  principally  of  four  NIR  absorbers  Hb-R,  Hb-02,  H2O,  and  lipids.  The 
concentrations  of  these  absorbers  (c)  were  quantified  by  solving  the  equation  p^  =  E  c  ,  where  E 
is  a  7  X  4  matrix  that  contains  the  molar  extinction  coefficients  of  the  four  chromophbres  at  the 
seven  source  wavelengths  ([23]-[25]). 
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For  each  measurement  we  report  four  hemoglobin  parameters:  [Hb-R],  [Hb-02],  total 
hemoglobin  concentration  (THC  =  [Hb-R]+[Hb-02]),  and  the  tissue  hemoglobin  saturation 
(St02=[Hb-02]  /  THC  X  100%),  where  the  characters  []  denote  concentration  (M  L’*).  Values  for 
water  and  lipid  content  are  reported  as  percentages.  The  water  percentage  is  the  concentration  of 
measured  tissue  water  divided  by  the  concentration  of  pure  water  (55.6  M).  The  lipid  percentage 
is  the  kilograms  per  liter  of  lipid  measured  relative  to  an  assumed  ‘pure’  lipid  density  of 
0.9  Kg  U'.  Thus,  the  water  and  lipid  percentages  we  report  are  relative  figures  of  merit 
compared  to  pure  solutions  of  the  substance,  and  are  not  strict  volume  or  mass  fractions. 

The  scattering  properties  of  the  tissue  also  yield  important  physiological  information.  NIR 
scattering  in  tissue  has  the  following  dependence:  =  A  where  ^4  is  a  constant,  X  is  the 
wavelength  (nm),  and  SP  is  the  scatter  power.  Scatter  power  is  related  to  the  scattering  center 
size  (d)  compared  to  the  optical  wavelength.  As  an  example,  scatter  power  is  4  in  the  case  of 
Rayleigh  scattering  (d«X)  and  is  ~1  for  large  Mie-like  scatterers  {d  ~X). 
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RESULTS 


A.  Sensitivity:  Pre-  versus  Post-Menopausal 

'» 

Both  Figure  2  and  Figure  3  present  a  typical  series  of  measurements  of  seven  absorption  and 
scattering  coefficients,  respectively.  The  points  represent  an  average  of  several  measurements  in 
the  center  of  the  left  upper  outer  quadrants  of  two  volunteers;  a  32-year  old  pre-menopausal 
woman  (squares)  and  a  54-year-old  post-menopausal  woman  (triangles).  Error  bars  show  the 
standard  deviation  of  repeated  measurements.  There  are  vast  absorption  and  scattering 
differences  between  pre-  and  post-menopausal  breast  tissue.  The  solid  lines  of  Figure  2 
represent  a  weighted-least-squares  fit  of  the  seven  absorption  coefficients  using  published 
extinction  coefficients  for  Hb-R,  Hb-02,  H2O,  and  lipids  [23-25].  Lines  between  the  rheasured 
points  have  been  interpolated.  The  solid  lines  of  Figure  3  represent  a  fit  of  the  scattering 
coefficients  to  the  equation;  Ps’  =  ^ 

These  optical  spectra  provide  unique  insight  into  breast  tissue  composition.  There  are  higher 
concentrations  of  hemoglobin  (i.e.,  both  Hb-R  and  Hb-02)  in  the  pre-menopausal  subject,  as 
evidenced  by  the  overall  higher  absorption  in  the  670  to  850  nm  range.  There  is  also  more  water 
relative  to  lipids  in  the  pre-menopausal  subject,  as  revealed  by  the  large  water  absorption  peak  at 
980  nm.  Recovery  of  the  absorption  spectra  allows  calculations  of  the  tissue  concentrations  of 
Hb-R,  Hb-02,  H2O,  and  lipids.  In  addition,  the  light  scattering  intensity  is  significantly  lower  in 
the  post-menopausal  breast.  Table  1  provides  a  summary  of  the  fitted  physiological  properties 
for  these  two  subjects  represented  in  both  Figure  2  and  Figure  3. 
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B.  Quantities  vs.  Age 


The  next  three  figures  present  FDPM-measured  quantities  as  a  function  of  volunteer  age. 
Figure  4  shows  the  THC,  Figure  5  shows  the  tissue  water  percent  (referenced  to  the 
concentration  of  pure  water),  and  Figure  6  shows  the  scatter  power.  Premenopausal  subjects 
(i.e.,  age  <50)  display  a  variety  of  values  in  all  three  plots.  This  spread  is  the  result  of  inter¬ 
subject  variations,  including,  but  not  limited  to,  menstrual  cycle  differences  and  gynecological 
age. 

There  is  a  general  decrease  with  age  in  THC,  water,  and  scatter  power  for  peri-  and  post¬ 
menopausal  volunteers  (i.e.,  >  50  years  old).  The  decrease  in  THC  correlates  well  with  previous 
histological  studies  showing  both  the  atrophy  of  well-vascularized  lobular  tissue  and  the  increase 
of  the  fat-to-collagen  ratio  after  menopause  [15].  This  is  consistent  with  compositional  analysis 
data  showing  lower  blood  and  water  content  for  fat  versus  glandular  tissue  [26].  Scatter  power 
also  decreases  with  age  after  50.  Compared  with  fat,  collagen  and  glandular  tissue  scatter  light 
with  a  higher  intensity  and  a  steeper  spectral  dependence  [27].  Thus,  smaller  scatter  powers  are 
expected  in  fatty  tissue.  Exact  knowledge  of  the  weighted  biological  contributions  to  optical 
tissue  scattering  is  a  matter  of  debate. 


C.  Correlation  Between  Absorption  and  Scattering  Information 
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Figure  7  displays  a  correlation  plot  of  the  measured  water  (open  squares)  and  lipid  (triangles) 
percents  vs.  scatter  power.  The  lipid  percent  is  the  mass  density  of  lipids  in  the  tissue  referenced 
to  the  mass  density  of  soybean  oil.  Note  that  neither  the  water  nor  lipid  percents  represent  a 
strict  mass  or  volume  fraction  for  the  tissue  (i.e.,  all  components  are  not  normalized  to  100%). 
Lines  represent  fits  of  water  (solid)  and  lipids  (broken)  as  functions  of  scatter  power. 

Tissue  scattering  correlates  well  with  the  parenchyma  composition  measured  via  absorption.  As 
scatter  power  increases,  the  lipid  percent  decreases  (/?^  = -0.841  ,P<  10'^)  and  water  percent 
increases  (7?^  =  +0.850,  p  <  10'^).  The  absorption  spectra  of  water  and  lipids  overlap  which  can 
lead  to  errors  in  measuring  these  chromophores.  However,  Figure  7  clearly  shows  that  both 
scattering  and  absorption  provide  complementary  information  regarding  breast  composition, 
indicating  that  such  errors  do  not  significantly  distort  our  findings. 


D.  Tissue  Metabolism  and  Menopause 

St02  is  the  percent  fraction  of  arterial  and  venous  Hb-02  that  contributes  to  the  tissue  THC,  and- 
should  not  be  confused  with  the  Sa02,  measured  via  pulse  oximetry  (which  is  the  arterial 
saturation).  Table  2  presents  the  measured  St02  in  all  subjects  with  a  body-mass-index  below  40. 
The  St02  of  the  pre-menopausal  volunteers  and  the  HRT  volunteers  both  averaged  approximately 
76%.  Their  means  were  not  statistically  significantly  different.  However,  the  St02  values  for  the 
post-menopausal  women  not  taking  HRT  were  significantly  higher,  averaging  81.9  ±8.8% 

{p  =0.062). 
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Decreases  in  St02  are  associated  with  increased  metabolism.  Tumors  for  example  have  been 
shown  to  have  lower  St02  values  than  non-tumor  tissue  [18].  The  post-menopausal  increase  in 
St02  relative  to  the  pre-menopausal/HRT  women  is  plausible  since  general  tissue  metabolism  in 
the  breast  should  decrease  after  menopause  as  the  tissue  begins  to  atrophy.  The  similarity 
between  the  St02  values  of  pre-menopausal  and  HRT  women  suggests  that  their  metabolisms  are 
similar,  reflecting  the  effects  of  both  endogenous  and  exogenous  hormones. 
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DISCUSSION 


Although  it  is  difficult  to  directly  validate  non-invasive  in  vivo  measurements,  these  initial 
results  indicate  that  the  sensitivity  of  NIR  spectroscopy  is  a  reasonable  refleetion  of  long-term 
hormone-controlled  breast  remodeling.  Additional  work  in  the  field  has  demonstrated  that  short¬ 
term  menstrual-eycle  changes  are  detectible  using  a  similar  NIR  teehnique  [28].  Our  results 
provide  indireet  validation  of  the  general  accuracy  of  NIR  breast  spectroscopy.  It  must  be 
stressed  that  this  quantitative  physiological  information  is  not  obtainable  by  any  other  non- 
invasive  radiologie  method. 

An  important  feature  of  NIR  methods  is  the  ability  to  characterize  quantitatively  the  breast  tissue 
of  women  regardless  of  age,  hormonal  status,  or  mammographie  density.  Inereased 
mammographic  density  contributes  to  a  22%  false  negative  rate  as  well  as  a  high  false  positive 
rate  (56.2%  cumulative  risk  after  ten  exams)  in  women  less  than  50  years  of  age  [29,  30].  A 
recent  study  found  that  routine  initial  mammography  was  not  clinically  advantageous  for  women 
less  than  35  years  of  age  [31].  Furthermore,  the  use  of  HRT  in  post-menopausal  women  is 
known  to  increase  mammographic  density  [32]  and  has  been  recently  shown  to  impede  the 
effieacy  of  mammographic  screening  [33,  34].  The  sensitivity  of  NIR  spectroscopy  to  known 
biological  processes  suggests  optical  methods  may  provide  important  information 
complementary  to  conventional  diagnostic  techniques,  particularly  in  the  ease  of 
radiographically  dense  breast  tissue. 
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V 


In  general,  pre-menopausal  breast  is  more  optically  attenuating  than  post-menopausal  breast, 
both  in  absorption  and  scattering.  It  is  important  to  note  that  this  trend,  particularly  in  scattering, 
is  also  true  for  mammographic  density.  Radiographically  dense  breast  tissue  is  due  to  differing 
amounts  of  fat,  collagen,  epithelium,  and  water  [35].  Although  x-rays  are  inherently  sensitive  to 
dense  electron  structures  in  tissue,  such  as  micro-calcifications,  NIR  photons  are  naturally 
sensitive  to  tissue  constituents.  Thus  NIR  spectroscopy  has  unique  potential  for  quantifying  the 
elements  of  breast  tissue  that  contribute  to  mammographic  density.  This  observation  may  be  of 
importance  in  screening,  since  it  may  identify  breast  tissue  at  physiological  risk  for  malignant 
transformation  and  can  be  performed  easily  in  all  women. 

In  addition  to  scattering  contrast,  we  report  age-dependent  variations  in  the  THC,  H2O,  and 
lipids,  as  well  as  changes  in  St02  for  pre-  and  post-menopausal  subjects.  Knowledge  of  the- 
“normal”  values  of  these  chromophores  will  play  a  role  in  evaluating  the  usefulness  of  optical 
methods  in  detecting  and  characterizing  lesions  in  the  breast.  Several  investigators  have  reported 
a  2-4  fold  THC  contrast  between  normal  and  tumor  structures.  In  vivo  tumor  St02  values  are 
also  typically  lower  than  normal  tissue.  [18,  6,  36].  Thus,  detailed  studies  of  normal  tissue  are 
essential  for  determining  the  sensitivity  required  of  optical  instrumentation  for  detecting  lesions 
in  women  of  varying  age  and  hormonal  status.  Furthermore,  as  baseline  levels  are  characterized, 
data  on  an  individual’s  absorption  and  scattering  variations  could  provide  important  insight  into 
disease  appearance  and  progression.  When  applied  to  patients  receiving  chemotherapy  and/or 
HRT,  this  information  could  also  be  used  to  generate  feedback  that  would  permit  customized 
treatment  planning  based  upon  individual  physiologic  response. 
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Table  1  -  Measured  physiological  properties  of  a  pre-menopausal  and  a  post-menopausal  breast 


[Hb-R] 

[Hb-02] 

St02 

THC 

[H20]‘‘ 

LIPID*’ 

SP 

(m) 

(laM) 

(%) 

(llM) 

(%) 

(%) 

PRE 

12.6±0.7 

27.9±2.7 

68.9  ±  1.3 

40.4  ±  2.7 

45.1  ±3.6 

29.8  ±5.6 

0.864  ±  0.068 

POST 

2.36±0.49 

12.2±1.6 

83.7±2.5 

14.4±1.9 

10.3±0.8 

67.5±3.7 

0.555±0.036 

^  Measured  concentration  relative  to  pure  water 

^  Measured  mass  density  relative  to  pure  lipid 

Hb-R:  reduced  hemoglobin 

Hb-Oa:  oxygenated  hemoglobin 

St02:  total  hemoglobin  saturation 

THC:  total  hemoglobin  concentration 

SP:  scatter  power 

PRE:  pre-menopausal  volunteer 

POST :  post-menopausal 
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Table  2  -  Comparison  of  tissue  hemoglobin  saturations 


S,02  (%) 

Observations  (#) 

PRE 

76.3  ±  6.4 

15 

HRT 

76.9  ±  9.2 

6 

POST 

81.9  ±8.8 

6 

PRE:  pre-menopausal 

POST:  post-menopausal 

HRT :  hormone  replacement  therapy 

StOi:  total  hemoglobin  saturation 
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Figure  1  -  Schematic  drawing  of  FDPM  instrument,  hand-held  probe,  and  measurement  map. 
The  components  of  the  instrument  are:  (1)  diode  lasers  (2)  avalanehe  photodiode  (3)  network 
analyzer  (4)  steady-state  current  source  (5)  bias  network  (6)  optieal  switeh.  Measurements  were 
performed  with  a  hand-held  probe  (close-up  on  bottom  of  Figure)  that  contained  the  APD  and  a 
source  optical  fiber.  The  measurement  location  was  in  the  center  of  the  left  upper  outer  quadrant 
as  indicated. 

Figure  2  -  Measured  absorption  spectra  for  32  year-old  pre-menopausal  (squares)  and  54  year 
old  post-menopausal  (triangles)  volunteers.  Points  represent  the  average  of  several 
measurements.  The  lines  represent  a  least-squares  fit  (extrapolated  to  all  wavelengths)  assuming 
the  breast  absorption  is  due  to  only  Hb-R,  Hb-02,  H2O  and  lipids.  Concentrations  of 
chromophores  from  this  measurement  are  listed  in  Table  1. 

Figure  3  -  Measured  seattering  spectra  for  same  volunteers  in  Figure  2.  The  lines  represent  a  fit 
to:  \is'  =A  where  S  is  the  seattering  (mm"'),  A  is  a  constant,  X  is  the  wavelength  (nm)  and 
I  SP I  is  the  scatter  power.  Results  of  the  fit  from  this  measurement  are  listed  in  Table  1 . 

Figure  4  -  Plot  of  the  total  hemoglobin  concentration  (micro-molar)  versus  volunteer  age 
(years).  Pre-menopausal  volunteers  (<50  years,  N=\5)  show  considerable  variation  with  age 
because  of  inter-subjeet  variations  such  as  menstrual  cycle  variations  and  overall  hormone 
production  differences.  Post-menopausal  (N=  10)  and  peri-menopausal  {N=  3)  volunteers  (>  50 
years)  show  a  marked  decrease  in  total  hemoglobin  concentration.  {N=  28  total  subjects.) 
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Figure  5  -  Plot  of  the  measured  tissue  water  percentage  (%)  versus  volunteer  age  (years).  Pre¬ 
menopausal  volunteers  show  considerable  variation  over  age  because  of  intra-subject  variations 
such  as  menstrual  eycle  variations  and  overall  hormone  production  differences.  Post¬ 
menopausal  and  peri-menopausal  volunteers  show  a  marked  decrease  in  water  eontent. 
Percentages  are  water  concentrations  normalized  to  the  concentration  of  pure  water  (N-  28  total 
subjects.) 

Figure  6  -  Plot  of  the  measured  optical  scatter  power  (arbitrary)  versus  volunteer  age  (years). 
Optical  scatter  power  refers  to  the  magnitude  of  the  slope  of  the  scattering  vs.  wavelength  curve. 
Pre-menopausal  volunteers  show  considerably  higher  scatter  power  than  the  post-menopausal 
volunteers.  This  general  pattern  is  the  same  for  mammographic  density.  {N  =  28  total  subjects.) 

Figure  7  -  Correlation  plot  between  water  percent  (open  squares)  and  lipid  percent  (triangles) 
vs.  scatter  power.  Percentages  are  referenced  to  the  concentration  of  pure  water  and  the  mass 
density  of  pure  fat.  Scatter  power  correlates  with  water  (/?^  = -f 0.850,  p<10'^)  and  anti¬ 
correlates  with  lipids  (R^  =  -0.841,  p  <  10'^).  (N=  28  total  subjects.) 
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Abstract 

Near  infrared  (NIR)  diffuse  optical  spectroscopy  (DOS)  and  diffuse  optical  imaging 
(DOI)  are  promising  methods  that  may  eventually  enhance  or  replace  existing 
technologies  for  breast  cancer  screening  and  diagnosis.  These  techniques  are  based  on 
highly  sensitive,  quantitative  measurements  of  optical  and  functional  contrast  between 
healthy  and  diseased  tissue.  In  this  study  we  examine  whether  changes  in  breast 
physiology  caused  by  exogenous  hormones,  aging,  and  fluctuations  during  the  menstrual 
cycle  result  in  significant  alterations  in  breast  tissue  optical  contrast.  A  non-invasive, 
quantitative  DOS  technique,  Frequency-Domain  Photon  Migration  (FDPM),  was 
employed.  Measurements  were  performed  on  fourteen  volunteer  subjects  (UC  Irvine  IRB 
#95-563)  using  a  hand-held,  deep-tissue  probe.  Optical  properties,  i.e.  intrinsic  tissue 
absorption  (pa)  and  reduced  scattering  (ps’)  parameters,  were  calculated  from  FDPM  data. 
Wavelength-dependent  absorption  (pa  at  674,  803,  849,  956nm)  was  used  to  determine 
tissue  concentration  of  oxyhemoglobin  (Hb02),  deoxy  hemoglobin  (Hb),  and  total 
hemoglobin  (Hbtot=  Hb02  +  Hb);  tissue  hemoglobin  oxygen  saturation  (St02  =  100  x 
(Hb02/  Hbtot))  and  bulk  water  content.  Results  show  significant  and  dramatic  differences 
in  optical  properties  between  menopausal  states.  Average  pre-menopausal  pa  and  p ’s 
values  at  each  wavelength  are  2.5  to  3-fold  higher  and  16-28%  greater,  respectively,  than 
absorption  and  scattering  for  post-menopausal  subjects.  Absorption  and  scattering 
properties  for  women  using  HRT  are  intermediate  between  pre-  and  post-  populations. 
Physiological  properties  show  differences  in  mean  total  hemoglobin  (7.0  pM,  1 1.8  pM 
and  19.2pM)  and  water  concentration  (10.9%,  15.3%  and  27.3%)  for  post-,  HRT,  and 
pre-  subjects,  respectively.  Because  of  their  unique,  quantitative  information  content. 
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diffuse  optical  methods  may  play  an  important  role  in  breast  diagnostics  and  improving 
our  understanding  of  breast  disease. 
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Introduction: 


X-ray  mammography  is  widely  used  for  breast  cancer  screening,  the  most 
common  form  of  cancer  in  women.  Because  of  uncertainties  associated  with  radiographic 
density,  mammography  has  up  to  a  22%  false  negative  rate  as  well  as  a  high  false  positive 
rate  (56.2%  cumulative  risk  after  ten  exams)  in  women  under  50  years  of  age  ’  .  A 
growing  body  of  evidence  supports  the  notion  that  mammography  cannot  accurately 
distinguish  between  benign  and  malignant  tumors  in  younger  women^.  A  recent  study 
found  that  routine  initial  mammography  was  not  clinically  advantageous  for  women 
under  35  years  of  age"*.  Furthermore  the  use  of  hormone  replacement  therapy  (HRT)  in 
post-menopausal  women  is  known  to  increase  mammographic  density^  and  has  been 
recently  shown  to  impede  the  efficacy  of  mammographic  screening®’  Techniques  such 
as  magnetic  resonance  imaging  (MRI)  and  ultrasound  (US)  are  used  only  as  secondary 
procedures  due  to  factors  as  high  cost  and  poor  specificity  (MRI)  or  low  sensitivity  (US). 

Currently,  invasive  procedures  such  as  fine  needle  aspiration  or  surgical  biopsy 
are  implemented  in  order  to  provide  a  definitive  diagnosis.  Given  the  sub-optimal 
perfonnance  of  x-ray  mammography  in  pre-  and  peri-menopausal  women,  the  majority  of 
invasive  follow-up  procedures  are  performed  on  normal  or  benign  tissue  that  present  no 
malignant  disease*.  As  a  result,  the  use  of  non-invasive,  near-infrared  optical  methods  as 
a  supplement  to  present  techniques  for  diagnosing  and  detecting  breast  cancer  has 
generated  considerable  interest. 

Optical  methods  are  advantageous  because  they  are  non-invasive,  quanitative, 
relatively  inexpensive,  do  not  require  compression,  and  pose  no  risk  of  ionizing  radiation. 
A  promising  NIR  optical  technique  currently  under  development  is  frequency  domain 
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photon  migration  (FDPM).  FDPM  methods  have  successfully  detected  the  presence  of 
small  palpable  breast  lesions  in  vivo  in  women  with  previously  diagnosed  breast 
abnormalities^'’'. 

FDPM  employs  intensity-modulated  NIR  light  to  quantitatively  characterize  tissue 
in  terms  of  its  optical  parameters,  i.e.  the  reduced  scattering  coefficient  (ps’)  and  the 
absorption  coefficient  (pa)-  The  concentration  of  significant  NIR  absorbers  (deoxy-  and 
oxy-hemoglobin,  water  and  adipose)  can  be  calculated  using  measured  pa  values  ’  . 
Multiple  light  scattering  in  breast  tissue  occurs  as  a  consequence  of  spatial  variations  in 
refractive  index  which,  in  turn,  are  influenced  by  cellular  and  extracellular  matrix 
density'"’’’^.  Both  pa  and  ps’  together  provide  an  understanding  of  changes  in  tissue 
cellularity,  metabolic  activity,  physiology  and  host  response  to  cancer.  Detection  of 
lesions  is  based  upon  functional  contrast  between  normal  and  diseased  tissue  in  the  same 
patient.  However  the  physiology  of  healthy  breast  tissue  is  complex,  influenced  by 
multiple  internal  and  external  factors  such  as  menstrual  cycle  phase,  menopausal  status, 
exogenous  hormones,  lactation  and  pregnancy.  Consequently,  in  order  to  establish  a 
better  basis  for  optical  detection  and  diagnosis  based  on  differential  functional  contrast, 
the  optical  properties  of  normal  breast  tissue  must  be  carefully  examined  and 
characterized. 

There  has  been  increasing  attention  paid  to  the  study  of  normal  breast  tissue 
optical  properties  and  changes  that  occur  with  age,  menopausal  status,  and  hormone 
replacement  therapy'^'’*.  During  the  reproductive  years,  the  breast  is  comprised  mostly 
of  glandular  tissue.  With  menopause,  there  is  a  progressive  atrophy  of  glandular  tissue, 
mitotic  activity  slows  down  and  vascular  requirements  decrease.  Some  of  these  changes 
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have  been  monitored  qualitatively  by  mammography  and  magnetie  resonanee  imaging"^’ 

2 1 

.  Involution  of  breast  glandular  tissue  can  be  impeded  by  the  use  of  exogenous 
hormones,  such  as  estrogen  and  progesterone  .  Presently  up  to  20%  of  U.S.  women  use 
hormone  replacement  therapy,  (HRT)  to  alleviate  symptoms  of  menopause.  However, 
long  term  use  has,  under  certain  conditions,  been  linked  with  increased  risk  for  breast 
cancer  .  Magnetic  resonance  imaging  studies  further  suggest  that  fluctuations  in 
endogenous  hormones  may  alter  menstrual  cycle  patterns  and  impact  breast  cancer  risk.^'^’ 
25, 26, 27  addition,  the  phase  of  the  menstrual  cycle  has  been  shown  to  influence 
mammographic  accuracy  and  survival  rates  following  tumor  resection^*’ 

In  this  paper  we  provide  initial  results  examining  differences  in  both  optical  and 
physiological  properties  of  healthy  breast  tissue  of  fourteen  women.  Our  results  show 
that  the  population  studied  can  be  well-differentiated  based  on  quantitative  near-infrared 
spectroscopic  measurements.  Optical  contrast  is  based  upon  physiological  changes 
accompanying  menstrual  cycle  variations,  exogenous  hormone  levels  and  menopausal 
status.  We  believe  this  study  is  particularly  important  because  standard-of-care  screening 
using  x-ray  mammography  has  diminished  efficacy  in  two  of  the  populations  we 
examine;  i.e.  pre-menopausal  and  HRT  subjects.  The  sensitivity  of  our  technique  to 
known  biological  processes  suggests  our  methods  may  provide  important  information 
complementary  to  conventional  diagnostic  techniques,  particularly  in  the  case  of 
radiographically  dense  breast  tissue.  Ultimately,  this  may  enhance  our  understanding  of 
pathophysiological  changes  that  accompany  malignant  transformation  and  provide  insight 
into  processes  associated  with  increased  disease  risk. 
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Materials  and  Methods 


a.  FDPM  instrument 

A  portable,  multi-wavelength,  high-bandwidth  frequeney-domain  photon 
migration  instrument  has  been  designed  and  optimized  for  clinical  optical  property 
studies  (figure  la).  The  instrument  employs  multiple  diode  lasers  (1)  to  provide  visible 
and  near-infrared  light  at  six  wavelengths  (674nm,  803nm,  849nm,  894nm,  947nm  and 
956nm  or  980nm).  The  FDPM  instrument  is  described  in  detail  in  Tromberg  et.  al., 

1997' '  and  Pham  et.  al.,  1999^°.  A  hand-held  probe  has  been  designed  to  house  an 
avalanche  photodiode  (APD)  that  records  the  diffuse  light  signals  after  propagation 
through  the  tissue  (2).  This  probe  has  multiple  grooves  on  the  casing  to  position  source 
optical  fibers  a  fixed  distance  from  the  APD  (Figure  lb).  A  1 00-pm-diameter  graded 
index  optical  fiber  positioned  25mm  from  the  APD  detector  delivers  the  diode  laser 
output  to  the  tissue  surface.  The  network  analyzer  (3)  measures  the  phase  and  amplitude 
of  the  electronic  signal.  A  DC  current  source  (4)  is  mixed  with  RF  power  provided  by  the 
network  analyzer  in  a  bias  network  (5).  This  distributes  power  to  each  laser  diode  and 
produces  amplitude  modulated  light.  An  optical  switch  (6)  delivers  light  serially  from 
each  diode  to  the  tissue.  The  optical  power  launched  into  the  subject  ranged  from  10- 
20mW  for  each  wavelength. 
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b.  Measurement  technique 

Subjects  were  measured  in  the  supine  position.  Measurements  were  made  in  reflectance 
geometry  by  placing  the  hand-held  probe  on  the  skin  surface  with  light  pressure.  The 
phase  shift  and  amplitude  of  the  photon  density  wave  was  measured  at  each  wavelength 
for  201  modulation  frequencies  between  50  and  lOOOMHz  at  a  source-detector  separation 
of  25  mm.  The  range  of  modulation  frequencies  was  swept  repetitively  so  that  each 
amplitude  and  phase  value  represents  up  to  twelve  measurements.  Four  to  six 
wavelengths  were  utilized  for  the  measurements. 

c.  Subjects 

A  total  of  fourteen  subjects  were  studied:  subjects  Pre  1  through  Pre  6  were  pre¬ 
menopausal  volunteers  18  to  33  years  of  age.  Post  1  through  Post  3  were  post¬ 
menopausal  women  age  57-67.  The  remaining  five  volunteers  were  post-menopausal 
women  51-60  years  of  age  using  hormone  replacement  therapy.  HRT  1  and  HRT  5  were 
on  combination  HRT  (progesterone  and  estrogen)  and  HRT  2  through  HRT  4  were  taking 
estrogen-only  HRT.  All  HRT  subjects  have  been  taking  exogenous  hormone  for  at  least 
2  years.  Experiments  were  conducted  in  adherence  to  UC  Irvine  IRB  approved  protocols 
95-563  and  99-2183.  The  volunteers  were  healthy  with  no  known  breast  diseases.  After 
providing  informed  consent,  volunteers  filled  out  a  brief  questionnaire  which  surveyed 
pertinent  medical  history. 

When  using  four  wavelengths,  the  time  required  for  a  single  measurement  was 
approximately  ninety  seconds.  A  total  of  twelve  sites,  six  per  breast,  for  each  volunteer 
was  measured  during  volunteer  visits.  The  measurement  map  is  shown  in  Figure  la. 
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Measurements  were  made  in  each  of  the  four  quadrants  placed  approximately  halfway 
from  the  center  of  the  breast  to  the  edge  of  the  breast  laterally.  Additional  measurements 
were  made  on  the  superior  areolar  border  and  the  axillary  tail.  Calibration  measurements 
were  made  approximately  every  ten  minutes  on  a  tissue  phantom  of  known  optical 
properties. 

A  separate  subject  (Pre  6)  was  studied  at  two  different  points  in  her  menstrual 
cycle.  One  measurement  was  performed  near  ovulation  (day  14)  of  the  subject’s  regular 
28  day  cycle,  and  another  measurement  performed  during  the  luteal  phase  (day  25)  of  the 
menstrual  cycle. 

d.  Model 

The  Pi  approximation  to  the  Boltzmann  transport  equation^'  was  used  to  analyze  the 
frequency-domain  data  and  is  valid  when  probing  homogeneous  turbid  media  at  source- 
detector  separations  greater  than  1  cm  and  modulation  frequencies  less  than  a  few  GHz. 
The  model  employs  an  extrapolated  boundary  condition  for  a  semi-infinite  geometry^^. 
The  solution  to  the  P  i  approximation  relates  the  frequency-dependent  phase  and 
amplitude  to  the  optical  absorption  (pa)  and  reduced  scattering  (ps’)  coefficients  .  To 
determine  the  optical  properties  from  a  given  set  of  frequency-dependent  data  a 
Marquardt-Levenburg  minimization  algorithm  was  used  to  simultaneously  fit  the 
amplitude  A(co),  and  phase  O(co)  by  minimizing  the  difference  between  the  measured 

values  and  those  predicted  by  the  Pi  approximation. 

Physiological  properties  were  calculated  from  the  determined  pa  values  assuming 
that  principal  chromophores  in  breast  tissue  for  the  NIR  wavelengths  used  are  deoxy- 
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hemoglobin,  oxy-hemoglobin  and  water.  The  concentrations  in  tissue  of  the  three 
components  were  calculated  using  pa  at  four  wavelengths:  674,  803,  849  and  956nm  or 
980nm.  The  contribution  of  each  component  to  the  total  pa  at  a  given  wavelength  is 
represented  by  the  equation  2.303  (£i  Cj)  =  //a,-,  where  £,  is  the  extinction  coefficient  in 

units  of  cm^/mole  of  a  particular  chromophore  at  a  given  wavelength  A,,  and  c,-  is  the 
concentration  of  the  chromophore.  The  total  absorption  at  a  given  wavelength  from  the 
tissue  chromophores  is: 


e  V)[Hb]  +  e\Hbo2)[Hb02]  +  e\H20)[H20]  =  pa\ 


(1) 


Where  the  brackets,  [  ],  indicate  concentration.  The  matrix  representation  of  the  four 


equations  is'^: 
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Given  that  there  are  more  wavelengths  than  principle  chromophores  we  have  4  equations 
and  3  unknowns  and  there  exists  no  general  solution  for  c,-.  The  chromophore 
concentrations  are  determined  using  a  least  squares  solution  to  Ec  =  pa-  In  matrix 
representation  the  chromophore  concentration  is  given  by:  c  =  (E  E)'  E  Pa  where  E 
and  E'*  denote  the  transpose  and  inverse  of  the  matrix  E  respectively. 
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e.  Calibration 


A  measurement  of  the  phase  shift  of  intensity  modulated  light  at  source  wavelength  X 
propagating  through  a  turbid  medium  at  a  given  source-detector  separation  can  be 
expressed  as  follows: 

^meas  ~  ^medium  ^instrument  (3) 

In  order  to  extract  from  at  a  single  source-detector  separation,  must 

be  explicitly  evaluated  (Eq.  (3)).  Similar  arguments  can  be  made  for  amplitude  data  (A  ); 
however  is  determined  by  evaluating  measured  and  instrument  response  ratios.  The 
instrument  response  was  calculated  by  making  a  measurement  on  a  homogenous  phantom 
with  known  optical  properties.  A  siloxane  phantom  cast  from  flexible  silicone  (RTV 
615/700,  GE  Corporation)  and  titanium  diode  in  a  400mL  mold  was  used  to  calibrate  the 
instrument  response  in  phase  and  amplitude  measurements.  Optical  properties  of  the 
phantom  were  determined  by  the  two-distance  technique  described  in  Fantini  et.  al, 
1994^^  Reference  measurements  are  made  on  the  phantom  at  the  given  source-detector 
separation  during  the  time  of  the  subject  measurement.  Calibrations  were  made  at 
approximately  ten-minute  intervals  to  account  for  amplitude  drift. 

Results  and  Discussion: 

The  results  of  raw  data  and  simultaneous  fits  to  the  Pi  approximation  for  Pre  1 
and  Post  1  are  shown  in  Figure  2.  Clear  differences  are  detected  between  the  two  subjects 
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in  both  amplitude  and  phase  data.  In  order  to  establish  the  range  of  optical  property 
values  characteristic  of  normal  menopausal  states,  we  examined  eight  volunteers  (5  pre¬ 
menopausal)  not  receiving  exogenous  hormone  therapy.  Figure  3  shows  a  scatter  plot  of 
Pa  versus  ps’  at  all  four  wavelengths  for  each  subject.  For  comparative  purposes  data 
from  the  right  upper  outer  quadrant  measurements  are  shown.  Glandular  tissue  is 
concentrated  in  the  upper  outer  quadrant  of  the  breast  and  consequently,  more  lesions 
appear  in  this  region  than  any  other^^. 

In  general,  post-menopausal  women  have  substantially  lower  pa  and  Ps’  values 
compared  to  pre-menopausal  women.  Average  pre-menopausal  pa  values  at  each 
wavelength  (0.0049  -0.016  mm’’)  are  2.5  to  3-fold  higher  than  post-menopausal 
absorption  (average  pa  =  0.00 16-0. 0064mm'').  Pre-menopausal  women  have  16-28% 
higher  scattering  values  at  each  wavelength  (average  ps‘  =  0.78-1. 1mm'')  than  post¬ 
menopausal  subjects  (average  ps‘  =  0.67-0.86  mm'').  Since  measurement  uncertainties 
are  less  than  5%  of  the  optical  property  value^*^  these  constitute  significant  and  dramatic 
differences  between  menopausal  states. 

A  deeper  understanding  of  the  underlying  physiological  reasons  for  optical 
property  contrast  can  be  obtained  by  examining  ps’  spectra  (i.e.  ps’  vs.  X )  for  three 
individuals.  Results  shown  in  Figure  4  highlight  the  impact  of  estrogen  on  breast 
structure.  The  29-year  old  pre-menopausal  subject  displays  ps’  values  that  are 
consistently  greater  than  the  67-year  old  post-menopausal  volunteer.  Interestingly,  ps’ 
values  for  the  52-year  old  post-menopausal  subject  receiving  estrogen-only  HRT  fall 
directly  between  pre-menopausal  and  post-menopausal  data.  This  may  be  due  to  the  fact 
that  estrogen  increases  the  rate  of  mitosis  within  ductal  tissue.  Thus,  optical  property 
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differences  between  HRT  and  post-menopausal  subjects  could  be  a  consequence  the 
impact  of  epithelial  tissue  proliferation.  Scattering  contrast  between  pre-  and  post¬ 
menopausal  subjects  is  probably  due  to  a  combined  effect  of  the  loss  of  glandular 
epithelium  as  well  as  extracellular  matrix  remodeling.  The  more  gradual  ps’  vs  A.  slope 
shown  for  the  post-menopausal  subjects  reflects  large  particle  scattering  consistent  with  a 
high  percentage  of  fatty  parenchyma  found  in  this  age  group.  In  younger  women,  the 
relatively  steep  ps’  vs.  X  slope  is  likely  to  be  influenced  by  the  presence  of  extracellular 
collagen  in  addition  to  cellular/epithelial  factors. 

These  concepts  are  further  substantiated  by  Figure  5,  which  illustrates  an  inverse 
correlation  between  ps’  and  age  for  subjects  over  50.  Whether  the  trend  indicating 
increased  scattering  for  HRT  women  is  due  to  exogenous  hormone  use  is  not  yet  clear. 
Pre-menopausal  ps’  values  appear  to  be  independent  of  age,  most  likely  due  to  the  timing 
of  the  measurements.  The  five  pre-menopausal  subjects  were  measured  at  different 
points  of  their  menstrual  cycle,  which  can  significantly  affect  optical  properties. 

Average  tissue  hemoglobin  concentrations  (oxy-,  deoxy-  and  total)  and  tissue 
water  concentration  are  shown  in  Figure  6.  The  concentration  of  blood  vessels  and  blood 
flow,  indicated  by  total  hemoglobin  concentration  (THC  =  [Hb]  +  [Hb02]),  is  greatest  in 
pre-menopausal  breast  due  to  the  high  vascular  demands  of  the  tissue.  Blood  vessel 
density  and  blood  flow  decrease  as  a  result  of  menopausal  involution  of  glandular  tissue. 
This  effect  can  be  detected  from  the  significant  total  hemoglobin  concentration  difference 
between  pre-  and  post-menopausal  subject  groups. 

Women  using  HRT  exhibit  hemoglobin  concentration  values  that  are  intermediate 
between  pre-  and  post-menopausal  subjects.  Thus,  FDPM  measurements  maybe 
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sensitive  to  subtle  HRT  effects,  including  previously  documented  elevations  in  blood 
flow,  fibroglandular  volume,  and  epithelial  cell  proliferation  The  general  decrease 
in  total  hemoglobin  concentration  with  age  after  50  (Figure  7)  suggests  it  is  difficult  to 
determine  from  limited  existing  data  whether  the  measured  HRT  hemoglobin  effect  is 
real  or  simply  a  consequence  of  age-related  optical  property  changes. 

Figure  6  data  show  that  the  average  water  percentage  for  pre-menopausal  women 
(29  ±  6%)  is  more  than  a  factor  of  two  greater  than  the  mean  tissue  water  concentration 
for  post-menopausal  women  (1 1  ±  2%).  These  differences  reflect  the  high  water  content 
of  epithelial  connective  tissue  compartments  in  pre-menopausal  tissue  while  post¬ 
menopausal  breast  is  dominated  by  low-water  content  adipose.  Water  concentration  in 

HRT  subjects  is  slightly  greater  than  the  post-menopausal  group  (15  ±  7%).  Differences 

* 

may  be  due  to  hormone-induced  accumulation  of  fluids,  a  commonly  occurring  side  effect 
from  HRT^^ 

Investigating  positional  variations  in  optical  properties  between  pre-menopausal 
(Pre  1)  and  post-menopausal  (Post  1)  subjects  reveals  a  higher  degree  of  variation  in 
scattering  for  post-menopausal  women  (approximately  8%  for  all  wavelengths)  vs.  pre¬ 
menopausal  women  (approximately  4%).  These  values  were  calculated  by  normalizing 
the  standard  deviation  of  the  values  from  measurements  in  each  of  the  four  quadrants 
(n=4)  to  the  mean  (data  not  shown).  This  difference  may  represent  the  non-uniform 
glandular  involution  of  breast  tissue  that  accompanies  menopause  and  results  in  palpable 
differences  in  density  in  the  breast  tissue  of  post-menopausal  women'^®.  The  positional 
variation  in  the  absorption  coefficient  is  greater  for  wavelengths  corresponding  to 
hemoglobin  absorption  (674-849nm)  for  the  post-menopausal  subjects  (23-33%  vs.  17- 
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20%  for  pre-).  However,  wavelengths  corresponding  to  fat  and  water  concentrations 
(894,  947  and  956nm)  show  a  higher  variability  in  pre-menopausal  subjects  (13-16%)  vs. 
only  7-8%  in  post. 

Figures  3-7  strongly  suggest  that  physiological  changes  due  to  hormonal 
fluctuations  that  occur  over  a  period  of  many  years  can  be  detected  and  quantified. 
Previous  studies'^'  show  these  effects  are  also  detectable  within  the  menstrual  cycle  of 
pre-menopausal  women.  In  order  to  test  FDPM  sensitivity  to  menstrual  cycle  variations, 
we  examined  a  pre-menopausal  subject  during  ovulation  and  before  the  onset  of  menses. 
Table  2  data  summarizing  these  results  show  that  pa  and  ps’  are  higher  for  the  latter  part 
of  the  cycle  (Day25)  than  in  mid-cycle  (Day  14).  Absorption  differences  correspond  to 
the  calculated  physiological  properties.  Our  results  show  a  48.3%  increase  in  hemoglobin 
and  28.1%  increase  in  water  concentration  during  the  luteal  phase,  changes  that  are 
consistent  with  the  physiological  effects  due  to  ovarian  hormone  fluctuations  during  the 
menstrual  cycle.  After  ovulation,  blood  flow  to  the  breast  can  increase  by  up  to  50%, 
there  is  an  increase  in  breast  volume,  and  parenchymal  water  content  changes  by  an 
average  of  25%  during  the  latter  half  cycle^^’  Changes  in  scattering  were  3-5%  and  not 
statistically  significant. 
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Conclusions 


A  hand-held  photon  migration  probe  has  been  developed  that  ean  deteet  significant 
differences  in  absorption  and  scattering  properties  occurring  with  changes  in  age  and 
menstrual  cycle  phase.  The  sensitivity  of  quantitative  NIR  spectroscopy  to  breast  biology 
is  unique  among  radiologic  methods.  Consequently,  optical  techniques  may  eventually 
provide  a  practical,  non-invasive  means  for  enhancing  the  accuracy  of  tumor  diagnostics 
in  pre-  and  peri -menopausal  subjects  and  for  increasing  our  understanding  breast 
pathophysiology. 
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Figure  Captions 


Figure  1 :  (a)  Schematic  drawing  of  FDPM  instrument,  hand-held  probe  and  measurement 
map  of  healthy  volunteers.  The  components  of  the  instrument  are:  (1)  diode  lasers  (2) 
avalanche  photodiode  (3)  network  analyzer  (4)  DC  current  source  (5)  bias  T  (6)  optical 
switch.  See  text  for  detailed  description..  The  breast  is  divided  into  four  quadrants:  upper 
outer  (a),  upper  inner  (b),  lower  outer  (c)  and  lower  inner  (d).  FDPM  measurements  are 
made  in  each  quadrant,  the  areolar  border  (e),  and  on  the  glandular  tail  (f)  that  extends 
into  the  axilla. 

Figure  2:  FDPM  measurements  of  (a)  phase  lag  and  (b)  amplitude  versus  modulation 
frequency  for  a  pre-menopausal  woman  (Pre  1)  and  a  post-menopausal  woman  (Post  1). 
Source-detector  separation  =  2.5  cm,  wavelength  =  849nm.  Solid  lines  represent  best 
diffusion  model  function  fits  to  phase  and  amplitude  data  for  each  subject.  Error  bars  are 
on  the  order  of  the  marker  size. 

Figure  3:  Absorption  coefficient  pa  versus  ps’  for  five  pre-menopausal  and  three  post¬ 
menopausal  subjects  at  all  wavelengths.  Values  are  calculated  from  best  diffusion  model 
fits  to  phase  and  amplitude  data 

Figure  4:  Reduced  scattering  coefficient,  ps’  versus  wavelength  for  three  subjects  of 
varied  hormonal  and  menopausal  status.  Some  values  at  980nm  represent  values 
extrapolated  from  AA,'®  linear  regression  analysis  of  the  data. 
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Figure  5:  gs’  at  674nm  versus  age  for  all  subjects. 


Figure  6:  Mean  hemoglobin  concentration  (gM,  oxy-,  deoxy-,  and  total)  and  water 
concentration  relative  to  pure  water  (%)  for  each  subject  group.  Values  are  determined 
from  wavelength-dependent  absorption  values  at  674nm,  803nm,  849nm  and  956nm  or 
980nm.  Error  bars  represent  the  normalized  standard  deviation  to  the  mean  for  5  pre¬ 
menopausal,  5  HRT  and  3  post-menopausal  subjects.  Confidence  values  are  80.9%  for 
Hb,  86.7%  for  HbOa,  86.3%  for  THC  and  98.0%  for  H2O. 

Figure  7;  Total  hemoglobin  concentration  versus  age  for  all  subjects. 

Table  1 :  Calculated  physiological  properties,  hemoglobin  concentration  (pM,  oxy-, 
deoxy-,  and  total)  blood  oxygen  saturation  (%)  and  water  concentration  (M)  for  a  pre¬ 
menopausal  subject  at  Day  14  and  Day  25  of  a  28  day  menstrual  cycle. 
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Figure  3: 


Absorption  vs.  Scattering 
for  Post  and  Pre  subjects 
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figure  7 
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Table  1:  Variation  in  physiological  properties  during  the  menstrual  cycle. 


Parameter 

Day  14 

Day  25 

Hb  (pM) 

3.07  ±  0.04 

4.72  ±  0.05 

HbO,  (pM) 

7.56  ±0.04 

11.04  ±0.06 

Total  Hb  (pM) 

10.63  ±0.05 

15.76  ±0.08 

0,  Sat,  (%) 

71.12±0.53 

70.06  ±  0.53 

Water  (%) 

16.49  ±0.11 

21.12±0.15 
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